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CUBIC SPLINE FUNCTION INTERPOLATION IN ATMOSPHERE 


MODELS FOR THE SOFTWARE DEVELOPMENT LABORATORY: 

FORMULATION AND DATA 

By James C . Kirkpatrick 
Lyndon B . Johnson Space Center 

SUMMARY 


Six standard reference atmosphere models for the years 1962 , 1963 (Patrick 
Air Force Base) , and 1966 (for July and January at latitudes 30° and 60° N) are 
presented in tabular form . The tabulation , which is presented in block data 
format , is adequate for the accurate representation of the atmospheric parameters 
of pressure , density , speed of sound , and coefficient of viscosity as functions of 
altitude. The range of tabulated altitudes extends from 0 to 205 kilometers . Inter- 
polation for the desired parameters is performed by using cubic spline functions . 
The recursive relations necessary to compute the cubic spline function coefficients 
are derived and implemented in subroutine form . Detailed FLOWGM flow charts 
and FORTRAN listings of the atmosphere package are presented in the appendix . 


INTRODUCTION 


This report presents a tabulation in block data format of selected altitude- 
correlated values of pressure , density , speed of sound , and coefficient of viscosity 
for each of six atmosphei*e models . These values were selected such that they 
adequately represent their respective functions for accurate interpolation by cubic 
spline functions throughout the range of tabulated altitudes . Three companion 
subprograms, which form the preprocessor (subroutines SDAT and SPLNl) and 
processor (subroutine ATMSFL) , are presented . These subprograms, together 
with the data element , compose the spline fit atmosphere package . In the preproc- 
essor , the data for the desired atmosphere model are selected and processed for 
proper order eind unit consistency before SPLNl is used to compute the spline 
function coefficients . In the processor , the data table is searched to establish the 
interval in which the desired altitude lies , and then the values of the previously 
mentioned parameters are computed by using the cubic spline function coefficients 
for the appropriate altitude interval . 

The cubic spline technique provides an effective , easy-to-use method for 
the unerring reproduction of a function from tabulated data . The method also 
ensures continuity in the function and in its first and second derivatives . After 






the cubic spline function coefficients ha^ been computed , they may be stored for 
repeated use as long as the tabulated data remain unchanged . 


SYMBOLS 





Cj(k) 

02(k) 

Cafk) 

dA 

ds 

dx 

d0 

£ 

El 

ef 

gh 

hk 

I 
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edges of the beam 
defined in equation (37) 
defined in equation (38) 
arbitrary constants of integration 
defined in equation (39) 
defined in equation (68) 
defined in equation (69) 
defined in equation (70) 

differential element of a cross-sectional area of a beam 

differential element of arc length 

differential element of the length of the beam 

differential element of the angle through which a cross-sectional 
surface of a beam is rotated under load 

defined in equation (49) 

modulus of elasticity of the beam 

flexural rigidity (can be taken as unity for any plane curve 
without loss of generality) 

original unstressed length of the fiber 

typical fiber of the beam 

defined in equation (1) 

area moment of inertia 

distance from the neutral surface to a deformed fiber of a 
loaded beam 
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M(x) 



y(k) 

y(x) 

y' 


y'(x) 

e 

P 

<P 


bending moment as a function of the distance x measured from 
the end of the beam 

bending moment at the kth section of the beam or an approxima- 
tion to the second derivative at the kth section 

load that acts upon the beam 

defined in equation (44) 

support reactions of the beam 

stress 

arc length of a curve 

defined in equation (46) 

shear force at the kth section of the beam 

defined in equation (24) 

defined in equation (23) 

slope of the deformation curve 

defined in equation (22) 

strain (change in length per unit of length) 

radius of curvature of the neutral surface 

angle that the tangent to a curve makes with the axis of 
abscissas 


SPLINE FIT ATMOSPHERE MODEL PACKAGE 


Model History 

The spline fit atmosphere package discussed in this report, was originally 
developed for the space vehicle dynamic simulator (SVDS) program in Novem- 
ber 1972. In February 1973, the package was incorporated with appreciable 
modification, into the Shuttle optimal abort program (SOAP) . The package 
presented here is from an early version of the SOAP . For implementation into the 
SOAP, the SVDS version was modified, at the expense of storage, to improve its 
performance and execution requirements . The spline fit atmosphere package has 
performed accurately in both programs . 
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Description of the Models 

The spline fit atmosphere package consists of six atmosphere models: 

1 . The 1962 standard reference atmosphere 

2 . The 1966 standard reference atmosphere for July at latitude 30® N 

3 . The 1966 standard reference atmosphere for January at latitude 30° N 

4 . The 1966 standard reference atmosphere for July at latitude 60° N 

5. The 1966 standard reference atmosphere for January at latitude 60° N 

6. The 1963 Patrick Air Force Base (AFB) reference atmosphere 

Each model consists of a tabulation of altitude-correlated values of the atmospheric 
parameters of pressure , density , speed of sound , and coefficient of viscosity . All 
six models are correlated to the same table of altitudes and tabulated in the fore- 
going sequence . All tabulated atmospheric parameter values are expressed in the 
meter-kilogram-second (MRS) system of units. Altitude values are given in 
meters , pressure values are given in millibars , density values are given in kilo- 
grams per cubic meter , speed of sound values are given in meters per second , 
and coefficient of viscosity values are given in kilograms per meter per second. 
Each atmospheric parameter is represented as a function of altitude by 123 data 
points . These data have been selected on best-fit considerations , with particular 
emphasis on those regions in which the tabulation shows that the function is dis- 
continuous in its first derivative . Each model is represented in the 0- to 
205-kilometer altitude range. The tabulation was compiled from references 1 to 3. 

Because the tabulation for the 1966 atmosphere extends only t' 118 kilo- 
meters, the last data point is repeated to 205 kilometers for model consistency. If 
the model altitude bound<.ries are exceeded, the package returns the last tabulated 
parameter values at the violated boundary . 


Assumptions and Approximations 

The spline fit atmosphere model package was developed to accommodate the 
needs and restrictions of the SVDS program . As a result , temperature , an impor- 
tant atmospheric parameter , was not included in the list of parameters modeled 
when the package was built . This parameter can be added to any model desired 
at the cost of tabulating the parameter values required . 

The tabulation was performed in the MRS system of units because this 
system was common to aL the references used. A common system of units is advan- 
tageous because a single set of conversion factors can be used to place all models 
in any desired working units . The British engineering system of units was used 
to convert the atmospheric data in this package (i.e. , pressure in pounds-force 
per square foot , density in slugs per cubic foot , speed of sound in feet per second , 
and coefficient of viscosity in pounds-mass seconds per square foot) . 
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Atmospheric parameter function values were selected at various altitudes . 
The increment of altitude was varied as a function of altitude as follows . 


Altitude range, km Altitude increment, m 


0 to 3 

250 

3 to 5 

500 

5 to 20 

1000 

20 to 160 

2000 

160 to 204 

4000 

204 to 205 

1000 


Additional data points were added to mark the location in which the atmospheric 
parameter functions undergo first-derivative discontinuities . These discontinuities 
occur only in the speed of sound and coefficient of viscosity functions , which are 
affected because of temperature variations- As a result, when temperature is 
added to the atmosphere models , the choice of altitudes should accommodate dis- 
continuities in the first derivative of the temperature functions . Because of the 
need to intersperse additional points , an ordering scheme is included in the pre- 
processor SDAT to place all altitude values and associated function values in 
ascending order of the independent variable . This order must be established 
before calling subroutine SPLNl to compute the cubic spline coefficients. 

The small altitude increment used at the extremes of the tabulation permitted 
the approximation of the first derivatives of the function at the initial and terminal 
boundaries by forming the difference quotient . This assumption of linearity at the 
extremes of the tabulated function is the only assumption made in the model . (The 
theoretical arguments on which the algorithm is based result in the additional 
assumption that the curvature of the function can be represented by its second 
derivative . The data sample selected is expected to comply with this assumption . ) 


DERIVATION OF THE RECURSIVE RELATIONS FOR COMPUTING 
THE COEFFICIENTS OF THE CUBIC SPLINE FUNCTIONS 


The mathematical spline is a formal analog of the draftsman's spline, a 
flexible beam device used to fair out a smooth curve between specified data points. 
The mechanical spline is anchored on the drawing board at each plotted data point . 
The curvature of the spline results from the bending moments applied to the beeim . 

In deriving the recursive relations of the spline algorithm , the elementary 
theory of elasticity is useful. Consider a simply supported beam (i.e. , ends 
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unrestrained) under the action of an applied load , however distributed . For an 
elementary analysis, the following assumptions are made: 

1 . Sections of the beam , which were originally plane , remain plane . 

2. The material of the beam is homogeneous and obeys Hooke's law (the 
deformation is proportional to the applied load) . 

3. The moduli of elasticy for tension and for compression are equal. 

4. The beam is initially straight and is of constant cross section. 

Figure 1 shows a uniform homogeneous beam under the action of the load P . 
The beam is simply supported and is in static equilibrium under the action of the 
support reactions Rj^ and Rg. The differential element of length dx, bounded by 

the edges a, b, c, and d, is bent or deformed under load into the exaggerated 
shape shown in figure 2 where c and d are deflected to c' and d' . The sec- 
tions ab and cd are shown rotated relative to each other through an angle d6 , 
but they remain undistorted in accordance with assumption 1 . The fibers in the 
region of ac are contracted , or shortened , in compression , whereas those in the 
region bd are elongated , or lengthened , in tension . Because the beam is uni- 
form and homogeneous in going from contraction to elongation , a region conceivably 
exists in which the fibers maintain their original unstressed length . This surface 
of unstressed fibers is referred to as the neutral surface . In the deformation of a 
typical fiber gh located at distance E from the neutral surface , the elongation hk 
is the arc of a circle of radius E , which subtends the angle d0 and is given by 
the relation 


hk = E de 


( 1 ) 


The strain e (the change in length per unit of length) is found by dividing the 
deformation hk by the original unstressed length of the fiber ef . 


_hk E de 

~“ir 


( 2 ) 


If the radius of curvature of the neutral surfp 2 e is p , the curved length of ef is 
equal to p d0 . The strain can then be expressed as 


E de _ E 

p d0 ~ p 


E = 


(3) 
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Because the material is assumed to be homogeneous and to obey Hooke's law, the 
stress (force per unit area) in the fiber gh is given by 


S =E = -B 
e p 


(4) 


where E is the modulus of elasticity defined as the ratio of the stress S to the 
strain e. Clearly, the sign of S varies with the sign of B: negative or compres- 
sive stress for B negative and positive or tensile stress for B positive . 

The beam in figure 1 is shown in static equilibrium . As a result , the alge- 
braic sum of all forces and moments acting on the beam must be zero . The bending 
moment M (x) must be balanced by the resisting moment . The resisting moment 
is the sum of all stresses S acting on each differential area element of cross sec- 
tion dA multiplied by its moment arm B measured from the neutral axis as shown 
in figure 3 . Hence , the bending moment is 

M(x) = J^fiCS dA) (5) 


Replacing S with EB/p from equation (4) gives 


M(x) =EI/p 


( 6 ) 


where I is the area moment of inertia . 


1 ^ M(x) 
p El 


(7) 


The reciprocal of p , the radius of curvature of the beam , is called the curvature of 
the curve . For plane cui^es , the curvature is defined as the rate at which the tan- 
gent to the curve turns compared with the description of arc . Therefore . if (p is 
the angle that the tangent to a curve makes with the axis of abscissas and s is the 
element of arc length , then the curvature is defined as 


1 _ dm 
p ds 


( 8 ) 
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The slope of the tangent to the curve is 


y' 


dx 


tan (p 


or 


<p = arctan y' 


( 10 ) 


and 


dtp = 




( 11 ) 


The differential element of arc length ds is defined as 


ds = 



( 12 ) 


Therefore , 


1 _ 

p “ ds 



( 13 ) 


and 


(i 



( 14 ) 


The assumption is usually made that the slope of the deformation curve y* is small 
2 

so that y' can be ignored compared to unity . When this assumption is made , the 
differential equation of the deformation curve can be written as 


Mt>:) = EIy"(x) 


( 15 ) 
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The assumption concerning the elastic behavior of a beam made at the beginning 
of this discussion can be realized by a mathematical function of one variable . The 
assumption that the slope of the curve is small and that its squa/e can be neglected 
can be realized by an arbitrary mathematical function only over relatively short 
interveds . Under this condition , the curvature can be approximated by the second 
derivative of the function . This point must be remembered when selecting a set of 
data points to represent a function for spline fit . Sufficient data points must be 
chosen in the appropriate concentration such that in the regions where the curva- 
ture is large , each section of the curve can be considered as a simply supported 
oeam (adequately supported so that the slope of the deformation curve ^s negligible). 

The following definitions are required for the discussion that follows. 

1 . The shear at a section of a beam is the algebraic sum of all external 
forces acting on one side of the section . 

2 . The bending moment at a section of a beam is the algebraic sum of the 
moments of all the external forces on one side of the section . 

Hence , if a mathematical function is represented in a region of definition by n data 
points, the region of definition will have been divided into n - 1 sections (fig. 4) . 
If each section of the curve between given data points is considered to be a beam 
deflected under load , then the bending moment in the kth section between the kth 
and (k + l)th data points will be given by the equatic ' 


M(x) = r dx + Mj^ 


where and are the shear and bending moments at the kth section . 
Therefore , 


Clearly , where x = x^ , 


M(x)=V^(x-x^). 


M (Xj^) = and where 


“k 


or 


(17) 


“kn=''k(Vi“’ht)*”k 


(18) 


Solving for g^ves 


„ “k.i - “k 

Vi ■ \ 



( 19 ) 



1 


N 


Substituting equation (19) into equation (17) gives 


M(x) = 


( Vi ■ M ^ 


= M. 


__ V ' J ' y V 




( 20 ) 


Substituting the bending moment M(x) by its equivalent from equation (15) (the 
product El is known as the flexural rigidity and can be taken as unity for any 
plane curve without loss of generality) in equation (20) gives 


(* - .. (Vi ■ *) 


M, 




( 21 ) 


where AXj^ = ■ Xj^ • Integrating twice gives 


(* ■ .. (’^+1 ■ ’‘f 




2 Ax, ^ ^1 


(22) 


(x - Xj^)^ - xf 

y(x) = M. ^ + M. ^ - 

^ k+1 6 AXj^ k 


+ C^x + C 2 


(23) 


where and C 2 are constants of integration . At x = x^, 
x = Xj^^j, y(Xj 5 +i) = y^+i- Substituting yj^ and yj^^j in equation (23) gives 


Ax, 


yk 


= M, 


+ + C 2 


and 


Ax, 


^k+i = “k+i 


+ C 


iVi 


+ c. 


(24) 


(25) 
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Subtracting equation (24) from equation (25) gives 




Solving for from equation (26) gives 


(“k*l ■ "k)**k 


Aktl **k»l^k\ / 

V“k ® /V 


^k VM 


Substituting for Cj^ from equation (27) in equation (25) and solving for gives 


*^2 = *k*l 


l^Vi “k^*k\ ^kn "kti^^kl 
6 6 ) 


Substituting the expressions for and C 2 from equations (27) and (28) in 
equations (22) and (23) gives the expression far the cubic spline function in the 


kth interval as 


(x - Xfe) (Xfe+i - x)' 

= “k»l ^ 6Ax^ ~ * “k ' 6 aj^ - 




11 



T 


} 


I 


and its first derivative is given by 


y’(x) =M, 


(* ~ *k) (*k+l ~ *) 


k+1 2Ax, 




Equations (21) , (29) , and (30) provide expressions for the values of the 
cubic spline function and its first and second derivatives in the kth interval . These 
expressions can be written for the cubic spline function and its first and second 
derivatives in the (k - l)th interval by replacing k with k-1. Therefore, 
in the (k - l)th interval. 


(* ■ \-i)‘ 


1a; 


, briL 

Vl 6Ax, , 


/^k-l “k-l^’'k-2^, 

te r— jK- 


y’(x)=Mj^ 


' Vl)' 


(*k - *)' 


k-1 2AJ 


/ yk _ ^k^^k-i\ /yk-i _ ^-i^vi\ 

V^*k-1 ® / V^’hc-l ® / 


(* ■ Vi) (*k " *) 

= "k - -a t ^ ^-^ “k-1 
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A cursory inspection of equations (29) and (31) and equations (21) and (33) 
will show that neither the cubic spline function nor its second derivative will 
experience discontinuity in going from the (k - l)th to the kth interval . This may 
be shown simply be replacing x with in these expressions . To ensure that 

the cubic spline functions will experience no discontinuity in the first derivative 
in going from the (k - l)th to the kth interval , replace x with x. inequations 
(30) and (32) and equate the results. Thus, 


or 


“«k 


"^^k-1./ ^k 
2 \Ax^_i 




l\ _/^k-l ®*k-l^*k-l\ 

r\^*k-i “ 6 I 




^*k/yk+i “kn"^*k' 




AXi. 


)-( 


Ax, 6 


) 


(34) 


'Vl 


Ay, _ 


Ax._, + Ax. Ax, 

r- ^k-1 " ^ ^k " -6^»*k.l - Ax,. AX 


k-1 


k-1 


(35) 


where Ay^^ = y^^^j - yj^ and Ay^^ j = yj^ - Equation (35) permits the 

writing of k equations for the solution of k + 2 unknown values of M . Clearly , 
two second derivative values must be supplied at any two of the tabulated data 
points. However, if the second derivatives are not known, first derivative values 
may be used through equation (30) . Further , if both the first and second deriva- 
tives of the function are known at any one of the tabulated points , the set of simul- 
taneous equations given by equation (35) can still be solved. 

Solving for in equation (35) gives 

2(Ax,^ t Ax^_,)M,^ .. 6 ^ 

.1- k-1 Ax^\ax^ Ax^_ 

' \“k * ‘'k“k-l * 'k 
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where 


2(Ax^ + 


'"k-1 

Ax, 


6 Pk 


Assume a solution of the form 




so that in the (k - l)th interval 


“k-i = Vi“k"Vi 


Substituting equation (41) in equation (36) gives 


"k+l “k**k ^ *’k('^k- A * ’^k-l)* *^k 
' A**>kVl)“k*Vk-l-"°k 


Comparing equations (40) and (42) gives 


Qj^ \^^^k^k-l 
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( 

I 


1 


and 


or 


U 


^ ■ ^k^k-1 °k 


(45) 


^k = -Qk(Vk-l"‘^k) 


(46) 


Substituting equations (37) , (38) , and (39) in equations (44) and (46) gives 


Ax, 


Ax, 


Qv = 


k 2(AXj^ ■^'^^k-l) ^’^k-l^k-l ^(^k+1 ■ ^-l) ^^k-l^k-1 


(47) 



(48) 


Equations (40) , (41) , (47) , and (48) are the recursive relations that will be 
used for the solution of the set of simultaneous linear equations given by equation 
(35) . Before this is done, equations (29) and (30) will be put in a more convenient 
form by using the following difference operators . 


'*k = ^ 

'‘k+i = 

Adding and subtracting Xj^ to the right-hand side of equation (50) gives 


(49) 

(50) 


<>k+i“*- 


’‘kti * ’'k 


*k = ““k ■ ^*k 


(51) 
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Substituting equations (49) and (50) in equations (29) and (30) gives 


\ Pk+1 ^k+l^*k 

= “k-i - “k 


( ^k “k'^’‘kV 

6 j\+l 


y.(x)=M -M 

y (,x; 2 Axj^ k 2AXj^ ® / \^*k 


^k ^k^*k^ 


M 5^-m 

K+1 2AXj^ ^ K 2 AXj^ AXj^ 6 


where AM^^ = - Mj^. 

From equations (49) , (50) , and (51) , when x = Xj^, then ® 
dj^^j = “AXj^; when x = then dj^ = AXj^ and dj^^^ = 0. Substituting these 

results in equation (53) gives, for x = x^^, 


'W) 


t 






Ax, 6 


r^C\n " 2m J 


and , for x = x. 


"(Vi) 


ktl Z 




'y, M, Ax. 
•'k k k 


Ay, Ax, 

AXj^ 6 V k k+1/ 
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The use of these recursive relations is described in the following six cases . 

Case 1: Solution of the set of simultaneous equations given by equation (35) 

when the second derivatives of the function are known at the initial and terminal 

tabulated boundaries . In this case, y”(x.)=M, and y"(x )=M . From 

11 n n 

equation (41) with k = 2, if Qj^ = 0, then By using these values for 

Qj^ and Uj^, equations (47) and (48) can be used to generate values of and 
, recursively , for all values of k ranging from 2 to n - 1 . After these values 
are available, then all values of Mj^ can be computed by using equation (40) 
written in the form of equation (41) 

“k = \ ««> 


and propagating the solution backwards with k ranging from n-1 to 2. As M 
and Mj are known, all n second-derivative values will be available . 

Case 2: Solution of the set of simultaneous equations given by equation (35) 

when the first derivatives of the function are known at the initial and terminal 

tabulated boundaries . In this case, y'(x,)=D, and y'(x )=D . Substitute 

11 n n 

k = 1 in equation (54) for from equation (40) and solve for to get 





If Qj = -1/2, then from equation (57) , 


U 


1 



(57) 


(58) 


By using = -1/2 euid the value obtained for from equation (58) , compute 
all values of Qj^ tnd Uj^ from equations (47) and (48) for values of k reinging 
from 2 to n-1. With the value of and thus computed, substitute 
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equation (41) in equation (55) for such that at x = x^, 


Ay . Ax , / \ 

r'(x \ = D = r — + 2_L fo M + U , + 2M 1 

V n) n Ax^_j 6 V^n-1 n-1 n-1 ni 


Soh in^ for from equation (59) gives 


.pziV^ , 

‘n-1 \ " "-I 


2+0 


With the value of M known , the solution can be propagated backwai d by using 

equation (41) with k ranging from n-1 to 1. Thus, M. will be computed as 
the last value . 

Case 3: Solution of the set of simultaneous equations given by equation (35) 
when the first and second derivatives of the function are known , respectively , at 
the initial and terminal tabulated boundaries. In this case, y* (x^) = and 

y”(x^) = M^. Proceed as in case 2 excluding the computation for M^, which is not 

needed . 

Case 4: Solution of the set of simultaneous equations given by equation (35) 

when the first and second derivatives of the function are known , respectively , at 

the terminal and initial tabulated boundaries . In this case, y'(x )=D and 

n n 

y"(Xj) = Mj^. Proceed as in case 1 and compute the value of as in case 2. 

Case 5: Solution of the set of simultaneous equations given by equation (35) 
when the first and second derivatives of the function are known at the initial 
tabulated boundary . In this case, y'(x,)=D+ and y"(x.)=M, . From equation 
(56) , if Qj = 1 , then 
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With the value of Qj^ = 1 and the value of given by equation (61) , compute all 

values of and from equations (47) and (48) for k = 2 to n - 1 . With 

these values of and Uj^ known , compute values for the second derivative from 

equation (40) by propagating the solution forward for k ranging from 1 to n - 1 . 
The last value thus computed is . 

Case 6: Solution of the set of simultaneous equations given by equation (35) 

when the first and second derivatives of the function are known at the terminal 

tabulated boundary . In this case , y' (x ) = D and y” (x ) = M . In equation 

nn nn^ 

(55), substitute Mj^ in equation (41) with k = n - 1 and solve for to get 


Ay . ^x 1 

y'fx ) = D = . ,M +U , + 2M 'I 

V n/ n ^*n-l ® \ n-1 n n-1 n) 

U , = - /—V M (? "Q A 
n-1 Ax i \ n Ax ,/ n\ n-1) 
n-1 \ n-1/ 


(62) 


If = 1. then from equation (62), 


U . = ^ Id - ^ -- -V 3M 

Vl \ " ^Vl/ " 


(63) 


Solving for and from equations (47) and (48) gives 


Vi = - 


Ax 


k-1 


U 


k- 1 Ax 


U 


k-1 


I AXk 

^ + 2(AXk ^ AXk_i) 

Ax, 

k'C 


Ax 


k-1 


Ax, 


Ql 


^ H*k+i 


^k-i 


- 6 


/^yk-1 ^yk\l 

K-1 "^"k/l 


(64) 


(65) 


With Q _i = 1 and U , given by equation (63), compute all values of Q, 

and from equations (64) and (65) for values of k ranging from n-1 

to 2 . With the values of and known , compute values for the second 

derivatives from equation (40) by propagating the solution backwards for k 
ranging from n-1 to 1. The last value thus computed will be 
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Clearly , the solutions for case 5 and case 6 could have been obtained from 
equations (54) and (55) . In the solution for case 5, can be obtained from 

equation (54) and y'(Xj^+j) f-om equation (55) . The solution can be propagated 

forward recursively with k ranging from 1 to n - 1 . In the solution for case 6 , 
can be obtained from equation (55) and y' (Xj^) from equation (54) . The solu- 
tion can be propagated backward recursively with k ranging from n - 1 to 1 . 


If equation (52) is written with replaced by its equivalent 


dj^^j “ ~ given by equation v51) , the result is 




M, 




6Ax, 




k-H 

Ax,. 




)■ 




“k'^’‘k 


) 


Ax. ' 6 / (‘^k “ ^*k) ■ ^1 




Ax. 


k 




\ ,2 ,3 

2 ^k 6Ax, \ 


( 66 ) 


Equation (66) is the equation for the cubic spline function coefficients in the kth 
interval . If there are n data points in the tabulation , there will be n - 1 intervals 
and n - 1 cubic spline functions to cover the range of tabulated data . 

For computation purposes , equation (67) is usually written as 
y(x) =yj^ + ^ 

where 


Cj(k) 

C2(k) 

C3(k) 







AM, 
k^ 

6A x, 
k 


( 68 ) 

(69) 

(70) 
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Equations (67) through (70) show that the coefficients are functions of the interval; 
however , the coefficients need to be computed only once . After the computation is 
made, the set pertaining to the interval in which the value of the independent 
variable lies can be brought in to compute the value of the dependent variable 
needed. Because there are three coefficients for each interval, if there are n data 
points, there will be n - 1 intervals and 3(n - 1) coefficients. 


CONCLUDING REMARKS AND RECOMMENDATIONS 


The atmosphere package for the 1962 and 1966 standard atmosphere models 
is recommended without reservation. When these models were built, it was possi- 
ble to conduct comparison checks with the layered versions of these models in 
250-meter increments. In both packages, these tests revealed errors that were 
properly corrected. The package for the 1963 Patrick Air Force Base reference 
standard model was not subjected to such strenuous tests because the layered 
version of this model was not available. However, with the other models so thor- 
oughly tested, the choice of altitudes was established as adequate, and the only 
requirement for the Patrick Air Force Base reference standard model was to ensure 
that the data were tabulated correctly . The laborious task of comparing (digit for 
digit) each tabulated number in the model printouts with the numbers in the 
reference source was performed as carefully as possible, and all errors found were 
corrected . As a result , the accuracy of this model is reasonably ensured . 

To test the accuracy of these models after they have been recoded for use in 
the Software Development Laboratory, the only requirement is to have the package 
reproduce the tabulated data. Because the interpolation is to be performed in a 
different set of units , the output must be converted to the units of the tabulated 
data for ease of comparison . The results should agree within an error of no greater 

_Q 

than 10 . Additional verification can be obtained by evaluating the parameters of 

the models at various altitudes and comparing the results with those obtained from 
established operational sources , such as the space vehicle dynamic simulator and 
the Shuttle optimal abort programs . These tests should be performed for all 
models . 

It is recommended that where the first and second derivatives are known at 
the same data point the cubic spline functions be used to interpolate data only in a 
small neighborhood of this point and not be extended for interpolation over the 
entire range of tabulated data . 


Lyndon B . Johnson Space Center 

National Aeronautics and Space Administration 
Houston, Texas, June 4, 1976 
986-16-00-00-72 
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APPENDIX 

PROGRAM FLOW CHARTS AND LISTINGS 


Definitions of the symbols used in the program flow charts and listings are 
as follows . 


SYMBOLS USED IN SUBROUTINE SDAT(IOP) 


lOP 


NF 

NP64 


NP34 


PRSLST 


an integer parameter used to designate the desired atmosphere based 
on the following set of values: 

lOP = 0 1962 standard atmosphere 

lOP = 1 1966 standard atmosphere for July at latitude 30° N 

lOP = 2 1966 standard atmosphere for January at latitude 30° N 

lOP = 3 1966 standard atmosphere for July at latitude 60° N 

lOP = 4 1966 standard atmosphere for January at latitude 60° N 

lOP = 5 1963 Patrick Air Force Base reference standard atmosphere 

an integer variable used to specify the number of data points used to 
describe each atmospheric parameter in each model (equal to 123) 

an integer variable used to specify the total number of data points used 
to describe all atmospheric parameters in each atmosphere model 
(equal to 6*NP or 6 X 123 = 738) 

an integer variable used to specify the total number of cubic spline 
coefficients for each atmospheric parameter (equal to 3*(NP - 1) 
or 3* (123 - 1) = 366) 

a column vector array dimensioned by NP64 or 738 containing all 
pressure data values (in millibars) for all atmosphere models 
tabulated 


RHOLST a column vector array dimensioned by NP64 or 738 containing all 

density data values (in kilograms per cubic meter) for all atmos- 
phere models tabulated 

SONLST a column vector array dimensioned by NP64 or 738 containing all speed 
of sound data \ les (in meters per second) for all atmosphere 
models tabulated 


A-1 



1 

1 ‘ 

- 1 1 1 ' - . 

1 ' { 

) 

1 

1 , ! : 

1 

VISLST 

a column vector array dimensioned by NP64 or 738 containing all coeffi- 
cient of viscosity data values (in kilograms per meter per second) 
for all atmosphere models tabulated 

ALTLST 

a column vector array dimensioned by NP or 123 containing the altitude 
data values (in meters) for all atmosphere models tabulated 

WORK 

a column vector array dimensioned by NP or 123 used as a working 
array for subroutine SPLNl 

D 

a column vector array dimensioned by 2 containing the values of the 
first (or second) derivative of the function at either or both of the 
tabulated terminal boundaries (only first- derivative values are used 
for the atmosphere package) 

ALTTAB 

a column vector array dimensioned by NP or 123 containing altitude 
values in the working units of the program (feet) 

PRSTAB 

a column vector array dimensioned by NP or 123 containing pressure 
data values for the selected atmosphere model in the working units 
of the program (pounds-foroe per square foot) 

RHOTAB 

a column vector array dimensioned by NP or 123 containing density 
data values for the selected atmosphere model in the working units 
of the program (slugs per cubic foot) 

SONTAB 

a column vector array dimensioned by NP or 123 containing speed 
of sound data values for the selected atmosphere model in the 
working units of the program (feet per second) 

VISTAB 

a column vector array dimensioned by NP or 123 containing the 
coefficient of viscosity data values for the selected atmosphere 
model in the working units of the program (pounds-mass 
seconds per square foot) 

CP 

a column vector array dimensioned by NP34 or 366 containing the 
cubic spline function coefficients for the interpolation for pressure 
in pounds-force per square foot 

CR 

a column vector array dimensioned by NP34 or 366 containing the 
cubic spline function coefficients for the interpolation for density 
in slugs per cubic foot 

CS 

a column vector array dimensioned by NP34 or 366 containing the 
cubic spline function coefficients for the interpolation for the speed 
of sound in feet per second 

A- 2 
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CV a column vector array dimensioned by NP34 or 366 containing the 

cubic spline function coefficients for the interpolation for the 
coefficient of viscosity in pounds-mass seconds per square foot 

J1 an integer variable used to designate the order of the derivative of 

the function supplied at the initial tabulated boundary (except when 
J1 = 4) . It can have the following values: 

J1 = 1 The first derivative of the function at the first tabulated point 
is supplied in D(l) . 

J1 = 2 The second derivative of the function at the first tabulated 
point is supplied in D(l) . 

J1 = 3 The first and second derivatives of the function at the first 
tabulated point are supplied in D(l) and D(2) , respectively. 

J1 = 4 The first and second derivatives of the function at the last 
tabulated point are supplied in D(l) and D(2) , respectively. 

JN an integer variable used to designate the order of the derivatives of 

the function supplied at the last tabulated boundary . It can have the 
following values: 

JN = 1 The first derivative of the function at the last tabulated point 
is supplied in D(2) . 

JN = 2 The second derivative of the function at the last tabulated 
point is supplied in D(2) . 

xs, ysl. Variables used for temporary storage in the ordering process for 
ys2, ys3,- altitude, pressure, density, speed of sound, and coefficient of 
ys4 viscosity, respectively 


SYMBOLS USED IN SUBROUTINE SPLNKN ,X,Y ,J1 ,JN ,D ,C ,W) 


N an integer variable used to specify the number of data points used to 

represent the function in the region of definition 

X a column vector array dimensioned by at least N in the calling element 

containing the values of the indepenaent variable in ascending order 

Y a column vector array dimensioned by at least N in the calling 

element containing the values of the dependent variable in sequence 
with the values in the X array 

J1 same as J1 defined in SDAT 
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JN 

D 

C 

W 

V 

IS 

FANS 


same as JN defined in SDAT 
same as D defined in SDAT 

a column vector array dimensioned by at least 3*(N-1) in the calling 
element containing the cubic spline function coefficient^ for the 
first- , second- , and third-degree terms for each interval . A portion 
of the array is also used to compute and store the values of one of the 
recursive functions used in the computation of the second derivatives 
8ind in the computation of the second- derivative values . 

a column vector array dimensioned by at least N in the calling element 
used as a working array to compute and store the values of one of 
the recursive functions used in the computation of the second deriva- 
tives and in the computation of the cubic spline function coefficients 


SYMBOLS USED IN SUBROUTINE ATMSPL(V.FANS) 


a real variable used to specify the altitude from the calling element 

containing the equations of motion at which values of the atmospheric 

parameters are desired 

an integer variable used to save the interval number in which the 

specified altitude lies 

a column vector array dimensioned by 8 containing the following: 

FANS(l) the interpolated value for pressure at V altitude in 
pounds-force per square foot 

FANS (2) the interpolated value for density at V altitude in slugs 
per cubic foot 

FANS (3) the interpolated value for speed of sound at V altitude in 
feet per second 

FANS (4) the interpolated value for coefficient of viscosity at V 
altitude in pounds-mass seconds per square foot 

FANS (5) the ratio c" the pressure at V altitude and the pressure 
at 0 altitude 

FANS (6) the ratio of the density at V altitude and the density at 0 
altitude 

FANS (7) the ratio of the speed of sound at V altitude and the speed 
of sound at 0 altitude 
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FANS (8) the ratio of the coefficient of viscosity at V altit'.de and the 
coefficient of viscosity at 0 altitude 

All other variables are defined in subroutine SDAT . 
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SUBROUTINE SOATCiOP) 

PARAMETER NP = 123. NP64 = 6«NP. NP34 » 3»NP - 3 
DAT A J I • JN| / I • I / 


c 

c 

c 

c 

c 

c 

c 

c 


lOP. 


.AN INTEGER VARIABLE RANGING FROM 0 TO 5 USED TO DESIGNATE 
THE ATMOSPHERE MODEL TO BE USED 

lOP = 0 - 19B2 standard ATMOSPHERE 

lOP = I - 19G6 standard for JULY AT 30 DEG. N. LATITUDE 

lOP = 2 - 1966 STANDARD FOR JAN. AT 30 DEG. N. LATITUDE 

lOP = 3 - 1966 standard FOR JULY AT GO DEG. N. LATITUDE 

lOP * 4 - 1966 standard fqr JAN. AT 60 DEG. N. LATITUDE 

lOP = 5 - 1963 PATRICK AFB REFERENCE STANDARD 

. ■■ ' 

SELECT DESIRED ATMOSPHERE INDEX ADJUSTQR 


{ lOPl = IOP*NP 


TRANSFER PERTINENT ATMOSPHERE DATA FROM STORAGE ARRAY TO WORKING 
ARRAY AND CONVERT TO DESIRED WORKING UNITS (BRITISH ENGINEERING: 





CdO 

1 = l.NP > 

- - - t>2 


|c 

TRANSFER 

altitude 

values list 

TO WORKING ARRAY 










[alt 

TABCn == 

ALTLSTC 1 1/0, 

.3048 1 






- ^ - 




TRANSFER 

PRESSURE 

values list 

TO WORKING ARRAY 









[ 

PRSTABCI] = 

PRSLSTCIOPl + n*2. 

11622/1.01325 









HI 

transfer 

DENSITY 

VALUES LIST 

• "T"" - “ ' ■' 

TO 'i^OR‘'.ING array 1 


E 


[RHOTABm = RHOLSTCIOPI + H »7 . 6474E-2/( 1 . 225*32. 1740491 

T 


TRANSFER SPEED OF SOUND VALUES LIST tq WORKING ARRAY 


SONTABCI) = SONLSTCIOPl ♦ 11/0.3048 
— 


CONT. ON PG 2 


SDAT 

PG .l QF. 
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ic TRANSTER COEFFICIENT OP VISCOSITY VALUES LIST TO trtJRKlMB ARRAY 



I CONTINUE I 

XS • ALTTAB(K) 

YSl « PRSTAB(K) 

YS2 « RHOTAB(Kl 
VS3 > S0NTA6(KI 
I VS4 « VlSTAB(K) 
ALTTAB(K) = ALTTAB(I) 
PRSTAB(K) * PRSTABU) 
RHOTABCK) :: RHOTAB(I) 


CONT. ON PC 3 




r 


Nv 


i 


SONTAB(K) = SONTAB(I) 
VISTAB(K) = \ 1STAB( 1 ) 
ALTTABU ) = XS 
PRSTABCIJ = YSl 
RHOTAB(l) = YS2 
SONTAB(I) = YS3 


2 - - - t> 




VISTABCI) 


VS4 


C APPROXIMATE THE FIRST DERlVATiVE OF THE FUNCTIQNCS) BY FORMING THE 
C FORWARD DIFFERENCE QUOTIENT OF THE INITIAL AND TERMINAL TABULATED 
C BOUNDARIES. THESE SLOPES ARE STORED IN OCI ) AND D(2). 

C RESPECTIVELY. 


E 


Join = (PRSTABC2) - PRSTABI in/(ALTTAB(2) - ALTTABUn 
D(21 = (PRSTABCNP) - PRSTA8CNP-m/CALTTABCNP) - A.LTTABTNP-l) ) 


COMPUTE CUBIC SPLINE INTERPOLATION COEFFICIENTS FOP PRESSURE (CP) 


(call SPLNKi ^.AlTTAB.PRSTAB.JI. JN.D.CP.WQrI^ 


D(l) = (RH0TAB(2) - RhOTABC I n/(ALTTAB(2) - ALTTAB(l)) 1 

0(2J = (RHOTAB(NP) - RH0TA6(NP-m/(ALTTAB(NP) - ALTTABCNP-ll) | 


fc COMPUTE CUBIC SPLINE INTERPOLATION COEFFICIENTS FOR DENSITY (CR)’ 
I CALL SPLNHNP. ALTTAB.RHOTAB. J1 . JN.Q.CR.WQRK) 


0(1) = (S0NTAB(2) - S0NTAB( 1))/(AlTTA8(2) - ALTTAB(D) 

0(2) = (SONTAB(NP) - SCNTAB(NP-1))/(ALTTAB(NP) - ALTTAB(NP- 1 ) ) 


C COMPUTE CUBIC SPLINE INTERPOLATION COEFFICIENTS FOR SPEED OF SOUND 
C (CS) 


(call SPLNKNP.ALTTAB.SONTAB.Jl.JNro.CS.WORK) 


cont. on pg a 


SOAT 

PG 3 OF 
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Pacp 





no 


1 




t 


( 


OCI) - (V1STAB(2) - VlSTAB(m/(ALTTAB(2) - ALTIABd)) 

0(2) s (VISTAB(NP) - V|STAB(NP-n)/(ALTTAB(NP) - ALTT ABCNP- 1 ) ) 


COMPUTE CUBIC SPLINE INTERPOLATION COEFFICIENTS FOR COEFFICIENT OF 
VISCOSITY (CV) 


ICALL yLNUNP»ALTTAB.VISTAB»Jl.JN.Q.CV.vbRM 


(return I 




SOAT 

eCLJ LIML 
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subroutine SPLNUN.K. Y. Jl. JN.D.C.W) 




CONT. ON P6 2 


SPLNl 

PG 1 or i 








non 
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no 
































c 


EMPLOY CUBIC interpolation tq COMPUTE OENSITy/SEA lEvEL RATIO 













i 



FAWS(2) = RHOTABCI) ♦ ((CRC3«n»DX ♦ CR( 3» I -I ) )»OX ♦ CRC3* I -2) 3»OX 
FANS(B) * FANS(2)/RHOTABCn 


^ ^ 


c employ cubic interpolation to compute speed of souno/sea level 

C RATIO 




FANS(3) » SONTAB(I) * CCCS(3*n»DX ♦ CS( 3» I -I ) )«OX ♦ CS(3*I-2) 3«OX 
FANS(7) s FANS(3)/S0NTABCll 


i - 

C EMPLOY CUBIC interpolation TO COMPUTE COEFFICIENT OF VISCCSITY/SEA 

C LEVEL RATIO 

i 


FANS(4) * VISTAB(I) ♦ (CCV(3«n*DX ♦ CV(3»I-in*DX ♦ CV(3» I -2) )«OX 
FANS(B) * FANSt41/VISTABm 



(return I 


|C INITIAL BOUNDARY VALUES - PUNCTIQ/SEA LEVEL RATIO \ 


td A2 


1 FANS( 1 ) 

1 

a PRSTABU) 1 


FANSC2) 

= RHOTABU) 

FANSC3) 

= SONTABm 

FANSC4) 

= VISTaB(I) 

FANSC5) 

= I.O 

FANS(6) 

= I.O 

FANS(7) 

= I.O 

FANSC8) 

» I.O 


(return I 


|c terminal boundary values - EUNCTION/SEA LEVEL RATIO 


- TT 

B 

3 


(FANSCn = PRSTABCNP) 

I "■■■ ■ ^ ■ ■ I 

FANSC2) = RHOTAB(NP) 

^ 

CONT. ON P6 4 ATMSPL 

efi-l- O F 4 


page is 

ijr POOR QlTAjr' 
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BLCCK DATA 
PARAMETER NP = 123 

PARAMETER MPII = t. NP|2 = IM. NP13 = 115. NPM * *tP 
PARAMETER NP2l=NPl UNP. NP22=NP12*NP. NP23=MP13*MP. NP2^=NPM*NP 
PARAMETER KP3t=NP2URP. KP32=NP22*NP, NP33=NP23*MP. NP34=NP24^NP. 
PARAMETER NP4l=MP3UKP» NP<2»NP32+MP. NP43=NP33*M». K'P44=NP34*NP 
PARAMETER N?5l=KP41*MP, NP52=t!P<2>NP, NP53=NP43*N?. NP54=NP44*NP 
PARAMETER NP61=M?3I*NP. NP62=NP54*S5. MPG3=NP54*9S. NP64=M?54*NP 












E 


ALTITUDE 

TABULATION 

IN METERS 1 












O 

> 

> 

tALTLSTtl). 

I=NP11.NP12)/ I 









ft 

1 

2 

3 

4 

5 

6J 




_4 




r 

O.OEO. 

2.5E2. 

5. 0E2. 

7.5E2. 

1.0E3. 

1.25E3. 









fc 

7 

8 

9 

10 

11 

. __‘2i 








f 

I.5E3. 

I.75E3. 

2.0E3. 

2.25E3. 

2.5E3. 

2.75E3. 1 








n 

13 

14 

15 

16 

17 

18 1 









|3.0E3. 

3* 31:3* 

4.0E3. 

4.5E3. 

5.0EE. 

6.0E3. 1 









E 

19 

20 

21 

22 

23 

24 1 









1 7. 0E3. 

8.0E3. 

9. 0E3. 

K0E4p 

1. IE4* 

1.2E4. 









G 

25 

26 

27 

28 

29 

5F1 









1 1.3E4. 

1.4E4. 

1.5E4. 

1.6E4* 

1.7E4. 

1.8E4. 1 








G 

31 

32 

33 

34 

3j 

36^ 









1 1.9E4. 

2.CE4. 

2.2E4. 

2.4E4. 

2.SE-?. 

2.8E4. 1 








Q 

: 37 

3Q 

39 

40 

41 

-‘Ll 




CONT. ON PG 2 PARTS 

PG \ CP 52 
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1.5E5. 1.52E‘i. 1.54E5. 


t03 104 105 


1.64ES. 1.68E5. 1.72E5. 


109 110 111 


1.S6E5. 1.S8E5. 1.60E5. 


106 107 108 


1.76E5. 1.80E5. 1.84ES. 


112 113 114 


I.88E5. 1.92ES. 1.96E5. 2. 0E5. 2. 04ES. 2.05E5/ 


DATA (ALTLSUn. I=NP13.NP14)/ 


115 116 117 118 119 i: 


121 122 123 


8.1E4. 8.3E4. 8.SE4/ 


. 1962 - STAUDARO ATMOSPHERE 

: PRESSURE VALUES TABULATED IN MILLIBARS 


DATA (PRSLST(1).1=NP11.NP12)/ 


1 2 3 4 5 6 


1.01325E3. 9.83576E2. 9.54612E2. 9.2G34SE2. 8.98762E2. 8.71G50E2. 


8 ' 9 


8.4558GE2. 8, 19988E2. 7.95014E2. 


. 13 14 15 


10 11 12 


7.70GG1E2. 7.4G917E2. 7.23771E2. 


16 17 18 


7.01211E2. 6.57803E2. G. 16G04E2. 5.77525E2. 5.40482E2. 4.7217GE2. 


ON P6 4 


















4.U052E2. 3.5GS16E2. 3.08007E2. 2.64399E2. 2.2S999E2. I.93994E2. 


1.6S79SE2. 1.41704E2. 1.2U18E2. 1.03528E2. 8.84971EI. 7.5G522EI. 


r^4S748El. 5.52930EI. 4.04749EU 2.97174EI. 2. 18837E1. 1.61619EI. 


1.19703E1. 8.890S3E0. 6.63412E0, 4.9P522E0. 3.77138E0. 2.87H3E0. 


2. 19967E0. 1.69496E0. 1.31340E0. 1.02296E0, 7.97790E-I.G.22283E-I. 


4.849i7E-l.3.7B572E-1.2.91373E-1.2.2460SE-l. 1.72457E-1. 1.31504E-1. 


9.94067E-2.7.44483E-2.5.52047E-2.4.05013E-2.2.9375SE-2.2. 1045E-2. 


I.4877E-2, 1.0386E-2. 7. IB91E-3. 4.9592E-3. 3.4313E-3, 2.374C"-3. 


I.6438E-3. 1.1449E-3. 8.0B83E-4. 5.74C3E-4. 4, 1377E-4. 3. 0075E-4, 



CONT. ON PG 5 
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2.2I23E-4. I.B500E-4. 1.2462E-4. 9.5225E-5. 7.3544E-5. 5.7623E-9. 


7S 80 81 82 83 84 


4.59I9E-5. 3.7137E-5. 3.0426E-5. 2.5217E-5. 7.1210E-5, 1.8139E-5. 


C 85 88 87 88 89 90 


1.5726E-5. 1.3791E-5. 1.2214E-S. 1.090SE-S. 9.0151E-U. 8.8882E-6. 


91 92 93 94 95 96 


8.0950E-6. 7.4104E-6. 6.8148E-6. 6.2931E-C. 5.8331E-6. 5.4252E-6. 


97 98 99 100 101 102 


5.0617E-6. 4.7345E-6. 4.4375E-G. 4.1G71E-6, 3.9202E-6. 3.6943E-6. 



103 

104 

105 

lOB 

107 

108 


3.2932E-6. 2.9475E-G. 2.6479E-B. 2.3851E-G. 2. 1536E-6. 1.9491E-8. 


, 109 110 111 112 113 114 


1.7680E-6. 1.B070E-6. 1.4330E-G, 1.3339E-G. 1.2179E-S. 1.1907E-B/ 


DATA (PRSLSTC I). 1=N?13.NPH)/ 


115 116 117 118 119 120 


3.31541E2. 2.09848E2. 3.46G86E;. 9. 033G7E-1.7. 04560E-1 .5.49540E-1 . 


121 122 123 


8.G204E-3, 5.9624E-3. 1.9758E-3/ 


. DENSITY VALUES TABULATED IN KILOGRAMS PER CUBIC METER 


ON P6 6 

















DMA (RHOLSTd). I=NPll.NPl2)/ 


2 3 4 5 


1.22SOEO. 1.1959E0. 1.16/3E0. 1.1392E0. 1.1117E0. 1.0846E0. 


8 9 10 11 12 


I.0581E0. 1.0321E0. 1.0056EO. 9.8151E-1. 9.S635E-1. S.3287E-1. 


13 14 15 16 17 18 


9.0925E-1. 8.6340E-1. 8.1935E-1. 7.7704E-1, 7.3643E-1. B.GOllE-1. 


19 20 21 22 23 24 


5.9002E-1. 5.257SE-1. 4.670BE-1, 4.1351E-1. 3.B480E-1. 3. 1194E-1. 



25 

26 

->7 

28 

29 

30 


2.GB60E-1, 2.2786E-1, 1.9475E-1. 1.6S47E-1, 1.4230E-1. 1.21B5E-1. 


31 32 33 34 35 3B 


1.0400E-1. 8.8910E-2. B.4510E-2. 4.B936E-2. 3.4257E-2. 2.507GE-2. 


1.8410E-2. 1.3555E-2. 9.8874E-3, 7.2579E-3, 5.3GGSE-3, 3.9957E-3. 


43 44 45 46 47 48 


2.9940E-3. 2.2589E-3. 1.7141E-3. 1.3167E-3. U 02S9E-3. 8.0097E-4. 


49 50 51 52 53 54 


B.3137E-4. 4.9762E-4. 3.9086E-4. 3. 05S2E-4. 2.3931E-4. 1.8837E-4 


ON PG 7 

















55 


56 


57 


58 


60 


t.47l3E-4. 1.1399E-4. G.7535E-5. 6.6593E-5. 5.0I51E-5. 3.73 6E-5. 


61 62 63 64 65 66 


2.750E-5. 1.999E-5. I.382E-5, 9.563E-6. 6.617E-6. 4.573E-6. 


67 68 69 70 71 72 


3.170E-6, 2.137E-6. 1.45SE-6. 1.008E-6. 7.0'14E-7, 4.974E-7. 


73 74 75 76 77 78 


3.493E-7, 2.492E-7. 1.804E-7. 1.323E-7, 9.82SE-8, 7. 153E-8. 


, 79 80 81 82 83 84 


5.321E-8. 4.035E-8. 3. ll2E-0. 2.436E-8. 1.844E-8. 1.434E-8. 


85 86 87 88 89 9 


1.I40E-8. 9.228E-9. 7.588E-9, 6.327E-9. 5.337E-9. 4.549E-9j_ 


, 91 92 93 94 95 96 


3.913E-9. 3.394E-9. 2.965E-9. 2.60DE-9. 2.307E-9. 2.053E-9. 


, 97 98 99 100 1 01 1 02 


1.836E-9, 1.G65E-9. 1.515E-9. 1.382c:-9. 1.264E-9. 1.159E-9. 


103 104 105 106 107 108 


9.970E-10. 8.624E-10. 7.532E-10. G.633E-10. 5.058E-10. 5. »09E-J^0, 




5 




C 109 

no 

111 

112 

113 

114 1 


ON P6 8 










4.608E>10. 4.115E-10. 3.632E-10. 3.318E-I0. 2.98GE-10. 2.9tOE-10/ 


DATA (RHOLSTU). 1=NP13.NP14)/ 


|C _ 

115 

116 

117 

118 

119 

120 


4.9576E-1. 3.3743E-1* 5.5006E-2. 1.1628E-3. 9.0690E-4, 7.1029E-4, 


. 121 122 123 


1.B62E-5, 1.15E-5. 3.810E-6/ 


. SPEED OF SOUND VALUES TABULATED IN MEiERS PER SECOND 


DATA (SONLST(I). I=NP11.NP12)/ 


2 3 4 5 6 


340.294. 339.333. 338.370. 337.403. 336.435. 335.463. 


89 10 11 12 


334.489. 333.511. 332.532. 331.549. 330.563. 329.575. 


D 

13 

14 

15 

16 

17 



19 20 21 22 23 


312.306. 308.105. 303.848. 299.532. 295.154. 


24 - 32 


9»295. 069. 


33 34 35 


ON PG 3 






















296.377. 297.720. 299.056. 300.386. 


37 38 39 40 41 42 


301.709. 303. 025. 306.489. 310.099. 313.665. 317.189. 


43 44 45 46 - 48 


320.672. 324.116. 327.521, 3*329.799, 


49 50 51 52 53 54 


327.911. 325.492. 323.058. 320.606. 317.630. 312.628. 


55 56 57 58 59 ' 60 


307.549. 302.307. 297.139. 291.800. 286.365. 280.03. 


, 61 62 - 114 


275.18. 53*269.44/ 


DATA (SONLSTC n. I=NP13.NP14)/ 


C 115 116 117 118 - 119 120 


305.984. 295.069. 297.049. 2*329.799. 329.114. 


C 


3*269.44/ 


C VISCOSITY COEFFICIENT 

C VALUES TABULATED IN KILOGRAMS PER METER PER SECOND 


DATA (VIbLSTC I). I=NPI 1.NP12)/ 


ON PG 10 































uu 


^ 61 62 - 114 


1.262E-S. S2*1.216E-S/ 


DATA (VlSLSTn).IsNPt3.K|Pl4)/ 


US tie 117 tie - its 120 


1.5099E-S. 1.4216E>5. 1.4376E-5. 2*1.7037E-5. 1.6381E-S. 


3»l.2l6E-5/ 


1966 STANDARD ATMOSPHERE FOR JULY AT 30-DEGREES NORTH LATITUDE 
PRESSIWE VALUES TABULATED IN MILLIBARS 


data (PRSLST(I), I=NP2l.NP22)/ 


2 3 4 5 S 


1.01350E3. 9.8S4E2. 9.579E2. 9.31CE2. 9.046E2. 8.788E2. 


Q.S36E2. 6.291E2. 8.05IE2. 7.817E2. 7.588E2. 7.385E2. 


13 14 IS 16 17 18 


7.148E2. 6.729E2. 6.331E2. 5.953E2. 5.593E2. 4.929E2. 


19 20 21 22 23 24 


4.331E2. 3.791E2. 3.307E2. 2.873E2. 2.48GE2, 2. 141E2. 


'25 26 27 28 29 30 


1.836E2. 1.S66E2. 1.329E2. 1.12SE2. 9.S25E1. 8.081E1. 


ON PG 17 













31 32 33 34 OT 36 


6.067EK 5.846EI. 4.258EI. 3.121EI, 2.300EI. t.705El. 


37 38 39 40 41 4^ 


I.270E1. 9.513E0. 7. 163E0. 5.425E0. 4.131E0. 3.162E0. 


C 43 44 45 . 46 47 48 


2.433E0, 1.881E0. 1.4G1E0, l.HOEO. 8.908E-1, 6.953E-1. 


49 50 51 52 53 ^ 


S.424E-1. 4.212E-I. 3.259E-1. 2.512E-1. :.924E-l, U462E-1. 


i — 


5 




IC ” 

56 

57 

58 

59 

~sT| 


1.102E-1. 8.225E-2. 6.083E-2. 4.45C~-2. 3.224E-2. 2.30GE-2. 


, 61 62 63 64 65 66 


1.629E-2. 1.135E-2. 7.792E-3. 5.301E-3. 3.607E-3. 2.455E-3. 


67 69 69 70 71 72^ 


1.671E-3. 1.142E-3. 7.877E-4. 5.488E-4. 3.859E-4. 2.737E-4. 


73 74 75 76 77 78 


1.358E-4. 1.427E-4. 1.065E-4. 8. 104E-5. 6.276E-5. 4.935E-5 


■Hi 

■■■■■■ 

■■■■•SHHHI 

IHHHHHHHH 

D 

79 

t 

o 

CD 

81 - 114 




3.955E-5. 3.240E-5. 

34*2.702E-5/ 


DATA CPRSLSn I ), 1=NP23.NP24)/ 


BN P6 13 





D 

115 

116 

117 

118 

119 

120 


3.542E2. 2.308E2. 3.643E1. 1.006EO. 7.874E-1. 6.146E-t. 


121 123 123 


9.423E-3. 6.427E-3. 2.025E-3/ 


. DENSITY VALUES T ADULATED IN KILOGRAMS PER CUBIC METER 


DATA (RHOLSTtI). I=NP2I.NP22)/ 


2 3 4 S 6 


I.ISSEO. 1.135EO. 1.1I2E0. 1.08SE0. 1.08GEO. 1.04IE0. 


8 9 10 11 12 


1.0I6E0. 9.923E-1. 9.6e6E-l. 9.452E-1. 9.222E-1. 8.997E-I. 


, 13 14 15 16 17 18 


8.776E-1. 8.348E-1. 7.937E-1, 7.540E-I. 7. 159E-1. 6.443E-1. 


19 20 21 22 23 24 


5.814E-1. 5.232E-1. 4.G94E-1. 4.199E-1. 3.742E-I. 3.324E-1. 


25 26 27 28 29 30 


2.941E-1. 2.592E-1. 2.275E-1. 1.929E-1. 1.617E-1. 1.357E-1. 


L 

31 

32 

' 33 

34 

35 

36 


I.142E-1. 9.618E-2. 6.870E-2. 4.944E-2. 3.579E-2. 2.607E-2. 


. 37 38 39 40 41 42 


ON P6 14 
















1.909E-2. l.40bE-2. 1.038E-2. 7.709E-3. 5.759E-3. 4.326E-3 


43 44 45 4G 47 48 


3.268E-3. 2.482E-3. 1.894E-3. 1.4G0E-3. 1.140E-3. 8.942E-4. 


49 50 51 52 53 54 


7.071E-4. 5.573E-4. 4.378E-4. 3.435E-4. 2.710E-4. 2.123E-4, 


55 56 57 58 59 60 


1.651E-4. 1.273E-4. 9.739E-5. 7.380E-5. 5.539E-5. 4. 113E-5. 


61 62 63 64 65 66 


3.020E-5. 2. 190E-5. 1.568E-5. 1.070E-5. 7.282E-6. 4.956E-6 


67 68 69 70 71 72 


3.374E-G. 2.25SE-6. 1.51GE-6, 1.030E-6. 7. 065E-7. 4.892E-7. 


73 74 75 76 77 78 


3.392E-7. 2.208E-7. 1.588E-7. 1.130E-7. 8.213E-8. 6. 087E-8. 


79 80 81-114 




1 4.300E-8. 3.259E-8. 

34*2.486E-8/ | 


DATA (RHOLSK n. 1=NP23.NP24)/ 


115 116 117 113 119 120 


i , 

4.95SE-I. 3.529E-1. 5.823E-2. 1.290E-3. l.OOSE-3. 7.955E-4. | 

' . ^ ■■ ■=; 

|C 121 122 123 


CONT. ON PG 15 


-32 









I.85GE-5. I.29QE-5, 4.089E-6/ 


SPEED OF SOUND VALUES TABULATED IN HEIERS PER SECOND 


DATA (SONLSTtn. |=NP2l.NP22)/ 


2 3 4 5 6 


349.9. 348.6. 347.3. 346.0, 344.7, 343.8, 


8 9 10 11 12 


342.9. 342.0. 341.1. 340.3, 339.4, 333.5, 


c 

13 

14 

15 

16 

17 

18 


337.7, 335.9, 334.2. 332.5, 330.7, 327.3. 


. 19 20 21 22 23 24 


322.9, 318.5, 314.0, 309.5, 304.9, 303.3. 


25 26 27 28 29 30 


295.6, 290.9. 286.0, 285.7, 287.2, 288.7, 


31 32 33 34 35 36 


290.2, 291.7, 294.6, 297.3, 299.9, 302.6, 


37 38 , 39 40 41 42 


305.2, 307.8, 310.8, 313.9, 316.9, 319.9, 


43 44 45 46 - 47 48 




— 5 — 



322.8. 

325.8. 

328.7, 

2t330.7, 

330.1, 


ON PG 16 































1 


2 


3 


5 


6 



ON PG 18 


M 















1.S07E-S. 1.467E-5. 1.427E-5. 1.30BE-5. 1.344E-5. 1.302E-5. 


, 61 62 63 64 - 67 


1.25SC-5. 1.21GE-5. 1.172E-5. 4tl.l68E-5. 


68 69 70 71 72 


1.190E-5, 1.218E-5. 1.245E-5. 1.273E-5. 1.299E-5. 


, 73 74 75 76 77 78 


t.335E>5. t.425E-5. t.S13E-5. 1.S99E-5. 1.681E-5. 1.762E-S. 


. 79 80 81 - 114 


_ 


I1.912E-S. 2.037E-5, 

34*2. 197E-5/ 


data (VISLS'.4).1=NP23.NP24)/ 


115 116 117 118 - 119 120 


1.593E-5. 1.483E-5. 1.429E-5. 2*1.71lE-5. 1.69SE-5 


; 121 122 - 123 


1.194E-5. 2*l.l6CE-5/ 


1986 STANDARD ATMOSPHERE FOR 
JANUARY AT 30-DEGREES NORTH LATITUDE 
PRESSURE VALUES TABULATED IN MILLIBARS 


DATA (PRSLSTC n. I=NP3l.NP32)/ 


2 3 4 5 6 


1.021E3. 9.912E2. 9.622E2. 9.340E2, 9.0S5E2. 8.V98E2. 



















8 9 10 11 12 


8.538E2. 8.285E2. 8.038E2. 7.798E2. 7.5B4E2. 7.335E2. 



13 

14 

15 

16 

17 

18 


7. 112E2. 6.682E2, 6.274E2. 5.88BE2, 5.517E2. 4.837E2. 


19 20 21 22 23 24 


4.22GE2. 3.679E2, 3.191E2. 2.757E2, 2.372E2. 2. 032E2. 


. 25 26 27 28 29 30 


1.734E2, 1.478E2. 1.257E2, 1.0S7E2. 9.040EU 7.G50EI. 


31 32 33 34 35 36 


6.479EI. 5.493EU 3.984E1. 2.907E1. 2.I34E1. 1.575E1. 


D 

37 

38 

39 

40 

41 



, 43 44 45 46 47 48 


2.194E0. 1.691E0, 1.310E0. 1.019E0. 7.941E-1. G.187E-1. 


4.808E-1. 3.723E-1. 2.872E-1. 2.207E-1. 1.GB7E-1, 1.281E-1. 


D 

55 

56 

‘ 57 

58 

59 

60 


9.G59E-2. 7.231E-2, 5.372E-2. 3.359E-2. 2.893E-2, 2.095E-2. 


61 62 63 64 65 66 


ON PG 2C 














I.503E-2. 1.068E-2. 7.538E-3. 5.323E-3, 3.7B0E-3. 2.G57E-3. 


. 67 68 69 70 71 72 


1.879E-3. 1.341E-3. 9.683E-4. 7.0S6E-4. 5.209E-4, 3.87GE-4. 


73 74 75 78 77 78 


2.910E-4. 2.210E-4. 1.699E-4. 1.322E-4. 1.039E-4. 8.239E-5. 


79 80 81 > 114 




1 6.606E-5. 5.;^5SE-5. 

34*4.394E-5/ | 


DATA (PRSLSTtn. I=NP33.NP34)/ 


115 116 117 118 119 120 


3.428E2. 2.196E2, 3.401EI. 8.997E-1. 7.010E-1. 5.457E-1. 


C 121 122 123 


8.970E-3, G.334E-3. 2.232E- 


OENSITY VALUES TABULATED IN KIl.OGRAilS PER CUBIC METER 


DATA (RHOLSTCI). I=NP31.NP32)/ 


2 3 4 5 6 


1.233E0. 1.200EO. 1.168E0. 1.137E0. 1.107E0. 1.078E0. 


8 9 10 11 12 


I.049E0. 1.021E0, 9.934E-1. 9.695E-1, 9.461E-1. 9.231E-1. 


3 14 15 16 


ON PG 21 




















73 lA 75 76 77 70 


4.322E-7. 3. 133E-7 , 2.305E-7. 1.71C[;-T. 1.29GE-7. 9.887E-8. 


79 80 81-114 


7.508E-8. 5.780E-8. 34M.509E-8/ 


DMA (RHOLSrm. I=NP33.NP34)/ 


115 116 117 118 119 120 


4.994E-1. 3.483E-1, 5.512E-2. 1.165E-3. 9. 073E-4. 7. 142E-4. 


121 122 123 


1.635E-5. 1.154E-5. 4.070E-6/ 


SPEED OF SOUND VALUES TABULATED IN METERS PER SECOND 


DATA (SONLST(I). I=NP31.NP32)/ 


2 3 4 5 6 


340.5. 340.0. 339.5. 339.1. 338.6. 338.1. 


8 9 10 11 12 


337.6. 337.1. 336.6. 335.6. 334.6. 333.5. 


: 13 

14 

15 

16 

17 



y 

1 332.5. 

330.5. 

328.5. 

326.4. 

324.4. 

320.3. 1 

'7— 

C 19 

20 

21 

22 

23 

2U 


316.2. 312.1. 307.9. 303.6. 299.3. 294.9. 


ON P6 23 PARTS 






























351.0. 360.3. 34*369.3/ 


DATA tSONLSTtn, I-NP33.NP34)/ 


116 117 118 - 119 







1 308.9. 

297.*. 

293.9. 

2*328.9. 

327.1. 1 


e 


3*2V7.2/ 


VISCOSITY COEFFICIENT 

VALUES TABULATED IN KILOGRAMS PER METEP PER EC^NO 


[data (VlSLSTt I). 1=NP31.NP32)/ 


2 3 4 5 6 


1.791E-5. 1.V87E-5. 1.783E-5. 1.779E-5. 1.V75E-5. 1.771E-5. 


8 9 10 11 12 


1.767E-5. 1.763E-5. 1.759E-5, 1.751E-5. 1.742E-5. 1.734E-5. 


13 14 15 16 17 1C 


1.726E-5. 1.709E-5. 1.693E-5. 1.67CE-5. 1.660E-5. 1.626E-5. 


Si 

19 20 21 22 23 24J 


1.593E-5, 1.559E-5. 1.525E-5, 1.491E-5. 1.456E-S. 1.420E-5. 


25 26 27 28 29 30 


1.405E-5. 1.391E-5. 1.37GE-f. 1.3G2E-5. 1.347E-5. 1.346E-5. 


31 32 22 34 35 36 


ON PG 25 PARIS 



















42 


I.359E-5. 1.373E-5. 1.401E-5. 1.423E-S. 1.445E-5. 1.466E-S 

37 38 39 40 41 



CONT. ON P6 26 PARTS 


ORIGINAL PAGE iS 
OF POOR QUALITY 


A-43 



on 


1366 STANDARD ATMOSPHERE PQR JULY AT 60-DEGREES NORTH LATITUDE 
PRESSURE VALUES TABULATED IN MILLIBARS 


DATA (PRSLSTin. I=NP4l,NP42)/ 


C I 2 3 4 S 6 


I.010E3. 3.804E2. 3.S16E2. 9.234E2. 8.960E2. 8.692E2. 


11 12 


8.432E2, 8. 177E2. 7.930E2. 7.68BE2. 7.453E2. 7.224E2. 



miimiiFvHiiim 


» 8 

9 

11 


13 

14 


15 

16 

It iFl 

■■■■■■■ 

■■■■■■■■ 

■■■■■ 




|7.000E2. 

6.571E2. 

6.164E2 

» 

5.779E2. 

5.414E2. 4.740E2. 




■HHUH 


BHHHBHl 


19 

20 


21 

22 

23 24 1 



■■■■■ 




I4.134E2. 

3.592E2. 

3.108E2 

» 

2.677E2. 

2.301E2. 1.978E2. 




HHHH 

HM 

HHliiHBlii 


C 25 

26 


27 

28 

29 30 1 


1.700E2. 1.4G1E2. 1.25BE2. 1.080E2. 9.285EI. 7.983E1. 


D 

31 

32 

33 * 

34 

35 

36 


6.864E1. 5.902EI. 4.3S4E1. 3.222EI. 2.398E1. 1.788E!. 


37 38 39 40 41 42 


1.338E1. 1.005EI. 7.593E0. 5.775E0. 4.421E0. 3.408E0. 


ON P6 27 PARTS 




















43 44 45 46 47 48 


2.640E0. 2.057EC, 1.607E0. I.259E0, 9.872E-1. 7.742E-1. 


49 50 51 52 53 54 


6.W2F-1, 4.751E-1, 3.704E-1. 2.877E-1. 2.218E-1. 1.694E-1. 


L ^5 56 57 58 59 60 


I.281E>1» 9.576E-2. 7.075E-2. 5. 16CE-2, 3.710E-2. 2.627E-2. 


61 62 63 64 65 66 


1.829E-2. 1.249E-2. 8.330E-3, 5.536E-3. 3.667E-3. 2.430E-3. 


87 68 69 7 0 71 72 


l.^lCE-3. 1.C75E-3. 7.297E-4. 5. 027E-4. 3.512E-4. 2.485E-4. 


73 • 74 75 76 77 78 


I.782E-4. 1.310E-4. 9.883E-5. 7.633E-5. 6. 004E-5. 4.801E-5. 



C 7 

9 8C 8 

1 > 114 





3.903E-5. 3.224E-5, 34*2.700E-5/ 


DATA (PRSLSTU ). 1SNP43.KP44)/ 


115 116 117 118 119 120 


3.343E2. 2. 133E2, 3.753EI. 1.115E0. 8.742E-1, 6.857E-I. 


C 121 122 123 


1.025E-2, 6.803E-3. 1.978E-3/ 


CONT. ON PG 28 






















7.K7E-4* 6.08SE-4. 4.8I4E-4. 3.81SE-4* 3. 049E-4. 2.416E-4. 


SS 56 57 58 59 60 


1.898E-4. 1.477E-4. 1.138E-4. 8.666E-5. 6.521E>5» 4.842E-5. 


61 62 63 64 65 66 


3.543E-S. 2.551E-5. I.800E-S* 1.192E-5. 7.898E-6* 5.233E-6. 



■■■■■■ 






c 

67 

68 

69 

70 

71 

72 


3.468E-6. 2.22SE-6. 1.454E-6. 9.640E-7. 6.49IE-7. 4.434E-7. 


73 74 75 76 77 78 


3.019E-7. 2.023E-7. 1.403E-7. 1.001E>7. 7.32<E>8. 5.46IE-8. 


79 80 81 • 114 


4.091E-6. 3.132E-8. 34*2.444E>8/ 


DATA £ RHOLST ( I) . I =NP43. NP44 )/ 


. 115 116 117 118 119 120 


4.S4IE-1. 3.301E-1. 5.807E-2. 1.40IE>3. 1.099E-3. 8.619E-4. 


121 122 123 


2. 150E-5. f.465F-5. 4.260E-6/ 


SPEED OF SOUND VALUES TABULATED IN METERS PER SECOND 


DATA CSONLSTI n. l=NP4I.NP42)/ 


2 3 


ON PG 30 





















61 62 


268.6. 261.9. 


68 69 


63 - 67 


5 * 255 . 0 . 


70 71 


259.9. 265.1. 270.2. 275.2. 280.1. 


73 74 75 76 77 78 


287.5. 301.1. 314.2. 326.7. 338.8. 350.8. 


m 


I 


365.5. 379.6. 34*393.2/ 


data CSONLSK I ). I=NP43.NP44)/ 


C US 116 - 117 118 - 120 


307.8., 2*300.8. 3*333.7. 


121 122 - 123 


258.3. 2*255.0/ 


VISCOSITY COEFFICIENT 

VALUES TABULATED IN KILCG'^AMS PER METER PER SECOND 


DATA (VISLST(I). 1=NP4I.NP42)/ 


2 - 3 4 5 6 


1.791E-5. 1.784E-5. 1.777E-5. 1.770E-5. 1.763E-5. 1.756E-5. 


8 9 to 11 12 


1.749E-5. 1.742E-5. 1.735E-5. 1.728E-5. I.721E-5. 1.715E-5. 


ON PG 32 






























1 


I 


1.14 IE-5. l.ie2E-S. 1.222E-5. 1.2G2E-S. 1.302E-S. 


73 74 75 76 77 78 


1.360E-5. 1.470E-5. 1.576E-5. 1.679E-5. I.777E-5. 1.87SE-5. 


: 79 

80 

81 - 114 

* — 

1 1.995E-5. 

2. lllE-5. 

34*2.222E-5/ | 


DATA tVISLST(I), I=NP43.NP44)/ 


115 116 - 117 118 - 120 


1.524E-5. 2»1.468E-S. 3*1.736E-5. 


121 122 - 123 


1.128E-5. 2»1.102E-5/ 


C ■ 


1966 STANDARD ATMOSPHERE FOR 
JANUARY AT 60-DEGREES NORTH LATITUDE 
PRESSURE VALUES TABULATED IN MILLIBARS 


DATA (PRSLSTCI). I=NP51.NP52)/ 


2 3 4 5 1 


1.01350E3. 9.804E2. 9.484E2. 9. 176E2. 8.878E2. 8.590E2. 


8.310E2. 8.038E2. 7.775E2, 7.519E2. 7.271E2. 7.031E2. 


13 14 15 16 17 18 


6.798E2. 6.352E2. 5.932E2. 5.534E2. 5. 158E2. 4.467E2. 


ON P6 34 











19 20 21 22 23 jfT 


3.853E2t 3.308E2. 2.823E2. 2.418E2. 2.067E2. 1.766E2. 


25 26 27 28 29 3 


l^0E2. 1.2dtE2. 1.103E2. 9.431E1. 8.058E1. 6.882E1> 


3» 32 33 34 35 36 


S.87SE1. 5.014E1. 3.647E1, 2.649EI. 1.922EI. t.338El. 


37 38 39 40 41 42 


1.020Elf 7.464E0> 5.479E0f 4. 041E0. 3. OOIEO* 2.243E0i 


3 44 45 46 


I.687E0. 1.277E0. 9.714E-1, 7.434E-1, 5.719E-1. 4.414E>1. 


49 50 51 52 53 54 


3.408E-1. 2.629E-1. 2. 02 IE- I . 1.540E-1. 1.184E-U 9. 050E-2. 


C 55 56 57 58 59 60" 


6.309E-2. 5.271E-2, 4.017E-2, 3.051E-2. 2.308E-2. 1.735E-2. 


61 62 63 64 65 66 


1.298E-2. 9.661E-3. 7.152E-3. 5.2G3E-3. 3.851E-3. 2.800E-3, 


. 67 68 69 70 ' 71 72 


2.022E-3. 1.4G0E-3. 1.06IE-3, 7.7G2E-4. 5.715E-4. 4.234E-4. 



73 

74 

75 

76 

77 

78 


ON P6 35 










f 



ON P6 36 

PARTS 


RS..35-QEL 


















2.422E-I. 2.071E-1. 1.770E-1, l.5l7E-lt 1.300E-1. I.U3E-1. 


31 32 33 34 35 36 


9.530E-2. e. 156E-2. 5.966E>2. 4.3b7t-2. 3. ir’»E-2, 2.275E-2. 


37 38 39 40 41 42 


U645E-2. t.l93£>2. 8.676E>3. 6.264E-3. 4.551E-3* 3.330E-3t 


43 44 45 46 47 48 


2.453E-3, 1.819E-3# 1.357E-3, 1.018E-3. 7.682E-4. 5.91 lE-4. 


O 50 51 52 53 54 


4.564E-4. 3.559E-4. 2.774E-4. 2. 150E-4. I.651E-4. 1.266E-4. 


55 56 57 58 59 60 


9.708E-5. 7.435E-5. 5.701E-5. 4.40OE-5. 3.392E-5. 2.599E-5. 


61 62 63 64 65 66 


^981E-5. 1.503E-5. 1.134E-5. 8.516E-6, 6.357E-6. 4.717E-6. 


67 68 69 70 71 72 


3,<80E-6. 2.467E-6. 1.V57E-6. 1.260E-6, 9.093E-7, 6.607E-7. 


73 74 75 76 77 78 


4.8r3£-7. 3.479E-7, 2.550E-7. 1.831E-7. 1.417E-7, 1.067E-7 





































1 


121 122 123 


298.5. 295.6. 286.7/ 


VISCOSITY COEFFICIENT 

VALUES TABULATED IN KILOGRAMS PER METER PER SECOND 


DATA fVISLSTin. I=NPS1.NPS2)/ 


. 1.639E-5. 1.642E-5. I.G44E-5. I.G47E-5. 1.643E-5. 


9 10 11 12 


1.639E-5. 1.634F-5. I.630E-5. 1.626E-5. 1.622E-5. I.6I8E-5. 


13 14 15 16 17 18 

1^ - — m: 

1.C14E-5. l.t. w-5^ 1.588E-5. 1.570E-5. 1.5b2E-b. I.5I6E-5. 


. 19 20 2 


1.480E-5. 1.443E-5. 7*1.424E-5. 


28 29 30 


1.421E-5. I.4I8E-5. 1.4I3E-5. 


. 31 32 33 34 35 36 


1.4I1E-5. 1.408E-5. 1.40IE-5. I.395E-5. 1.396E-5. 1.407E-5. 


C 

37 

38 

' . 39 

40 

41 

42 


1.418E-5. 1.429E-5. 1.440E-5. I.4GGE-5. I.492E-5. 1.519E-5. 


PG 40 
































i 

( 

• 

i 

t 

i 

j 1 

i 

1 

1 

i ' 

, f 

i i 


OAU (PRSLSTCl). |sNP6I.NP62)/ 


2 3 4 5 


1 . 0I70147EO. 9. 8829373E42. 9. 602265 1E«2. 9. 3280864E«2. 9. 0603418E« 2. 


6 7 8 9 10 


8.79895d8E^2.8.S4*JB573E^2.8.2949430E«2.8.052U68E»2.7.8152728E«2. 


11 12 13 14 IS 


7.5843002E«2.7.3590840E«2.7. 1395065E«^2.6.7167869E«^2.6.3iS1745E«2. 


16 17 18 19 20 


5. 933705 0E«2. 5. 57 1 4348E«^2. 4. 900S9 1 2£ ^2. 4 . 29879S9E«2. 3. 753204 0E«2. 


21 22 23 24 25 


3.2649869E«2. 2. 8277555E^2. 2. 4373U4E*2. 2. 090928tE^2. 1 .7861 068E«^2. 


26 27 28 29 30 


1 .5I99026E«2. 1 . 2892856E^2. 1 . C91 184 1E^2. 9. 22S2842E« 1 . 7.8097365E« 1 . 


31 32 33 34 35 


6. 62B0092E ♦ 1 . 5 . 63 1 5652E ♦ 1 . 4 . 0899 1 9 1 E ♦ 1 . 2 . 99 1 8759E ♦ 1 . 2 . 2 038 1 59E ♦ 1 . 


36 37 38 39 40 


1.6327363E»l. 1.2146273E>l.9.0905080E + C.e.84299l4E*0.5. 1807184E* 


41 42 43 44 45 


3.9447995E*0. 3. 0209180E*0. 2. 32624 1 lE-^O, 1 . 8004513E+0. 1 . 3994781E*0. 


ON P6 42 



































CONT. CN PG -14 PARTS 


A-61 








1.0S37666E«Q. I. 02eO2;E*0. t. 003S670E*0.9.790260l£»l«9.:> 4^»-* 


9.3122447E- I • 9. 0800345E> I • 8. 852568 IE-l»8.4t22243E-l»7.99l5662E- 


7.S9046^7E- » ♦ 2C34275E- 1 » 6. 4983435E- 1 * 5. 8535 153E- 1 . 5. 265 1817E-1 7 




4 .7749382~E- i . 4 > 2255460E- 1 « 3. 7638429E- 1 . 3. 3302 1 20E-I » 2. 92322 18E-I 




7.«i432637E-1.2.1920326E-i.l.C7l7665E-l.U5845601E-l7l73239218E-l. 


i.l096236E-l.9,3193799£-2.6»6193250E-2.4.747889eE-2.3.4382489E-2» 


f2.5119029E ■2>l,8334b60E>2.l.3457797E-2.9.930102eE-3.7.3654170E-3»' 


I S.4934199E-3.4. I220206£ ~3.3. U34715E-3.2.3684559E-3. U815l546E^ 


1 .40157eeE'-3. 1 . 0965534E-3. 8. 6526723E-4. 6. 825322 tE-4. 5. 37566B4E-4 


rp 

4.2227457E-4»3.3048920E-4»2»5745233E-4» 1 ♦9944483E-4» l«5352539E-4 


CONT. ON PG 45 


PARTS 






























































. 61 62 63 


2.7964961E2. 2.7570998E2. 2.7160535E2. 


64 - 67 


4»2.6944122E2. 






[C 

68 

69 

70 1 


■■■■■■ 




2.7387919E2. 2.7Q24640E2. 2.82S46I1E2. 


71 72 73 74 75 


2.8678t36E2. 2.9095497E2. 7.9778102E2. 3. 0445407E2. 3. 1098395E2. 


76 77 78 79 80 


3.1737952E2. 3.23G4874E2. 3.3d 83625E2. 3. 47ui)S71E2. 3.5897208E2. 


. 81 82 83 84 85 


3.6S3978SE2. 3.80704S1E2. 4.0126162E2. 4. 2081570£2» 4.3950065E2. 


. 86 87 88 89 90 


4.5742298E2. 4.74SS910E2. 4.9131 021E2. 5.0740585E2, 5.2300637E2. 


91 92 93 94 95 


5.3315483E2. 5. 528884 1E2. 5.B723942E2. 5.8123S20E2. 5.9490377E2/ 


data (SCNLST(I). I=NP63.NP64)/ 


96 97 98 99 100 


6.0826430E2. G.2l337b2E2, 6.309648BE2. G.4047440E2. 6.4979159E2, 


101 102 103 104 105 


B.590032GE2. G.6809798E2. B.8001215E2. G.9173079E2. 7.0153775E2 


ON P6 49 

































I 


i 


I 1.SS09244E-S. r.Sl0e096E-S. 1.4707842E-S. 1.433b282E-S. I.4008886E-S. I 
' T- - ■* 

|C 26 27 28 29 30 


I 1.3745403E-5. 1.3jS 8590E-S. 1.34S79S8E-S. 1.3447579E-S. t.3S85386E-5. I 

' ■ ' • — I— * 


E 


31 


32 


33 


34 


3S 


I 1.3724675E-1>. 1.3867977E-5. 1.4H5944E-S. 1 . 438£'j70E-5, 1.4587102E-5, I 

I , — , ] 

fc 36 37 38 3? 4 FI 

■ ^ , 1 

I I.4748972E-5. 1.4982730E-5. 1.522353SE-5. 1.5474 16C:-5, 1.5734080E-5. I 
• - 


[L 


41 


42 


43 


44 


45 


\ I.5999574E-5, 1.62G377Cc-5. 1.6S16764E-5. 1.6745577E-S. 1.69342C6E-5. | 


E 


46 


47 


48 


43 


50 




t 1.7063446E-5. 1.7034883E-5. 1.6870136E>5. 1.6660460E-5. 1.6414072E-5 
^ 


51 


52 


53 


54 


Hj 


[ 1.6138668E-5. 1.5841407E-5. 1.5528854E-5. 1.5206665E-5. 1.48B0579E-5 


56 


57 


58 


59 


iT] 


I 1.4554073E-5. 1.4230412E-5. 1.391I371E-S. 1.3597275E-5. 1.3286009E - 37~| 

^ 


E 


61 


62 


63 


64 - 67 


I 1.2976813E-5>1.2662082E-5tl.2335131E-5. 


4»1.2163172E-5. | 


E 


68 


69 


70 i 


1.2516126E-5. U286^SME-5. 1 . 3200747E-5. | 


CONT. P6 51 


PARTS 

Pf i 50 OF 52 


A-68 


f 





I 


I . 3548638E-S. I . 3884397E-5. I . 44351 OlE-5. 1 . 4975097E-5. ! . 5504842E-5. 


1.6024767E-5. 1.6535285E-5. 1.7529832E-5. 1.8490756E-5. 1.9421230E>5. 


2. 0323355E-5. 2.11991 97E-5. 2. 2873253E-5. 2. 4474 1 25E-5. 2. 59942 16E-5. 


2.7448117E-5.2.8842997E-5.3.0184890E-5.3. 147892eE-5,3.2729499£-5. 



3.394C400E-5.3.51H925E-5.3.6255?5SE-5.3.7388033E-5.3.844739iic-5/ 

, r.i: , 

DATA (VISLSK 1). 1=NP63.NP64)/ 


3.9502050E-5.4.0531788E-5.4. 1288666E-5.4.2033057E-5.4.2765514E-5. 


4.3486565E-5.4.4195706E-5.4.5127375E-5.4.60.404S5E-5.4.6803460E-5. 


4.7423143E-5.4.80352S3E-5.4.8640170E-5.4.9238014E-5.4.9745032E-5. 


5.01B3831E~5.5.0579335E-5.5. 0931617E-5.5. I 094193E-5. 1.5708688E-5. 


CONT. ON P6 52 


A-69 
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A-71 



I 




001^4 


C 

lOP . 2 • 1944 iTANOAKu ^09 JA||. At )Q 0C4. N. U»T!TuOC 

00IU4 


C 

1 

00104 

24* 

c 

lOP • 2 - 194* ST4N0AKU P09 JU^T At 40 dC4* M# i^AMTuOC 

0010* 

i7^ 

c 

* 

00104 

24» 

c 

top • • 1944 STAHOaRO 909 JAH. AT ^0 N« l.AT|T>jOC 

00 1 u* 

2»» 

c 


00104 

20* 

c 

lOP • 4 • 1942 PaTkIcA a 9» ACr^RrHCC STANOAHn 

OOlU* 

21* 

c 


ooto* 

22* 

c 


00104 

22* 

c 


00104 


c 

SCLCcT OCSIRCO aTmoSP>«£9C |wOC1 AdJU4T09 

CO 10* 

24* 

c 


0011 1 

24* 


10*1 • IOP*KP 

001 1 1 

22* 

c 


00 M 1 

24* 

c 

TRANSft« PCRTIncaT ATHOS^MCRC Ca’a 9R0H ST09a4C a99aT TO RQPaIng 

001 1 1 

29* 

c 

A99AT and COSWtR? TO OCStULO •0*A|NC UM|T$ n|9ltISH CN4tNCC9|N4l 

00 Ml 

90* 

c 


00 M< 

91* 


DO 1 1 « IfNP 

00 m2 

92* 

c 


00M2 

92* 

c 

TAANSrt* ALTITUoC VALUES tIST To A09K1NC a99«T 

001 i2 


c 


COM4 

94* 


AtTTASlM • AlTlSTM 1/0*2094 

ooms 

94* 

c 


00M4 

92# 

c 

TRAKlSrCR pRCSSUPC Vi^ocs tlST To f09K|N4 ARRAY 

00M4 

94* 

c 


00M4 

99* 


PRSTA4M» • P9S;.STM0RI * I) *2 • M 422/ 1 *0 1 22S 

00 M 4 

40* 

c 


00 M 4 

SI* 

c 

TRANSFi^ U£NS|TT VaCoCS UST TO RQf^lMNS *RR*T 

00M4 

42* 

c 


00 M^ 

42* 


RHOTaOMI • 9H0LSTM09I * M *T • 4979C»2/ M • 224* Jf* | 79099 > 

00 M ^ 

49* 

c 


00 M ^ 

44* 

c 

TRANSFtR S^eCO or -‘yRO VALOCS ttS7 Tq r0RK:N4 aRRAT 

00M7 

44* 

c 


00120 

42* 


SQHTAbMI • SONLSTMoRt * ll/0*209t 

00120 

44* 

c 


00120 

59* 

c 

transfer coefficient or viscosity vauucs list To rorain® array 

00120 

40* 

c 


00121 

41 * 


1 VISTABMI • YjSLSMInPl ♦ I l • 1 • 2029 / M • 749 9 • 22 • I 79 j » 

00(22 

42* 


H • NP • 1 

00122 

42* 

c 


00122 

494 

c 

RCA9RAN4E UaTa in ACcORoANCE RITH aSCCNqInN VALUES OF ThE 

00122 

44* 

c 

inoepenoenT Variable - altitude 

00122 

44* 

c 


00l2<« 

42* 


00 200 1 • 1 ,M 

00122 

44* 


A • ! 

00130 

49* 


IJ ■ 1 ♦ 1 

00121 

20* 


00 200 0 • tJ*NP 

0C12H 

2| » 


IF ( ALT r A8( A ) .OT* AlTTABU M A • J 

Q0124 

T2* 


200 CONTINUE 

OOMO 

22* 


A$ ■ ALTTaIMAI 

00 M 1 

29- 


rsi • PRSTABf 1 

00M2 

24* 


Y$2 • RHOTAB'AI 

00IH2 

2#* 


TSJ • SONTAB AI 

OOM** 

22* 


TS-* • V ISTA3 lA I 

001 '*4 

24* 


AlMAbI9> • ALTTf.BlM 

00 1 ' 4 

29* 


PRSr ABCL 1 ■ P9S' ABl 1 1 

0QM2 

42« 


KMJTABIAI * RHUMHMI 

001 

8 1 • 


40 IT AJ tK 1 * so* f A. i 1 1 


A-72 



OOlSI 

• 2* 


VISTASIR) • V1ST4BII 9 

QOIVl 

• 2* 


AUTTA4I I 1 • AS 

001^^ 



^AsrAtitl • rsi 

OOlSi 

• (• 


fiMOTAatn • TS2 

ooiss 

• 4« 


SQfVTABlIl ■ 

Q0IS4 

47* 

200 

VlSTASdl • TS4 

00IS4 

• •• 

C 


001^* 

• 4» 

c 


Q0IS4 

♦0* 

c 


00114 


c 

APF40I|MAT( The riAST 0 C*Ma7I»C Of ^h£ PUnc 7 1 0«« 1 9 1 ®T rOtMiNC TMc 

OOIS4 

♦ 24 

c 

r0B44A0 0|PFC«(NCt •UOTtCMT Of TMf lM|TtAL 4^0 tCO»ti«AL TaBULATCO 

001(4 

4J« 

c 

SOUmOAMICs* Th£SC Stores A«e STOACO IM Dll) ANO 0(21 « 

001(4 

44« 

c 

«CSPCC7I9CLT* 

001(4 

♦ (• 

c 


00l(( 

((• 

c 


00! (4 

4 7* 

c 


00140 



0111 • IPPSTltin * aBSTAOI 1 1 |/|AL^7A4lf 1 • ALTTAenn 

00141 



0121 • (PASTAtlAPl - P4$TA®tAP-| l|/IAtfTAtlMPI - AL T 7 AB fWA* 1 1 1 

0014k 

100* 

c 


00141 

101* 

c 

COAPuTl COBIC SPUINt INTCRPOLATiOn COtfFlClENTS FQn PMCSSuAC ICPI 

^0 i 4 I 

1024 

c 


00142 

102* 


CALL SPLNi (nP* alM AB«^^S7 A9 • J1 t JN , D t CP • «0Rk 1 

00142 

1044 

c 


00142 

IOC 

c 

* 

00142 

10(4 


0111 • CAhOTABIZI » BHOTABI 1 n/|ALTTAB(21 * ALTrABlUl 

0014*1 

I07« 


0(21 • (AhOTAB(nP| • AH07;.B(NP*|9|/(ALTrAB(NP} • ALT T AB | 1 1 1 

00144 

I0«4 

c 


00144 

Illf4 

c 

COAPUTC CMBIC SPLINE 1 NTEPPOL A T | On COcmriENTS ro* OLNSITY |CA| 

00144 

1104 

c 


0014( 

111* 


CALL SPLNI INP,ALTTAB«^H07ab. 2I • JN , 0 tCA tBOAK 1 

00I4S 

1124 

c 


OOI4( 

1124 

c 


00144 

1144 


0111 • IS0NTaB(2I • SONTAbC 1 11/(Al7TaB( 21 • ALTrABCIIl 

00147 

11(4 


0(2) • (SoNTAB(nPI - S0NTaB(NF* 1 1 |/(AtrTA0(NP) - AlTTAB |NP«| 1 1 

00147 

1144 

c 


00147 

117* 

c 

CONPuTE CutIC SPLINC I NTCAPoLA 7 | On COEfriCIEMTS FOB SFEtO OF SOu^O 

031(7 

1144 

c 

ICSl 

001(7 

119* 

c 


00173 

120* 


CALL SPCNt ( nP«AUTTAB.S0n7A4,3I « jN ,0 >CS •BOBA I 

COI 70 

121* 

c 


00170 

122* 

c 


00171 

I22« 


0(11 • (V1STAB(2> * VlSTABt 1 1 1 /(Al7Tab 121 • ALtTABMll 

00172 

:i4« 


0(21 • (VtSTABlNPl - V|$TAB(NP*1 ) |/( AlTTABInPI - AL T T A0 ( Np* C 1 1 

00172 

12(4 

c 


00172 

12«« 

c 

COMPUTE Cubic spline intebpolaT|On cOcfficients fob coefficient of 

00172 

1274 

c 

FISCOSITT (CVl 

00172 

1244 

c 


00172 

129* 


Call SPLNI (NP,ALTTAB»Vls7A0tUl • JN,0»CV»«0BA ) 

00172 

1204 

c 


00172 

1214 

c 


00174 

1224 


•'ETuBM 

00 » 

1224 


end 


iNU of 
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m r%m ii f 

IIM rMTt^s f tl k|*KW <*• *lK«tCt ktrit 

*■!• t0*i II AT 10|it|lt 


lM«t0T|«C AtMf^ C«T«T »tUr «t0ttl 


9T00AM iltfVI C0Oill) Mtff’l •ATalt) •M«<| t ALaKK C«»«90t<| 9O0900 
€•««•« •LViMi 
9««J 9^k»A* M«0>y 


| 9 |. 9 C«« VAaCi 


miW 


$T 9 aaAC AtSIAVMClir |A^ 0 C>| TT 0 C, tCkATlTt LSCATtt«t MNC| 


•••I •M0t« I0A. 

•001 0 M 0001 AkTtAI 
0000 0 000001 01 

000 J 0 e 00 i«« ■H 0 rA| 


1001 0C00» 100k 

0000 0 0OII«V 
0000 I 000001 t 

000 J « 000 »«l » 0 «Ta« 


0001 000000 took 

000 i ■ 00 irt 0 c« 

0000 000000 IMJ 99 

0009 • O00t»« «t>r*0 


0001 «00l>l 000k 

0009 0 1^.000 CS 
0000 I 0000 19 


0001 00OI0I 00lk 
0009 • 009 tAl Cf 
0009 0 O00IV1 900rtt 


00101 

1« 



1 


00101 

!• 

C 




00109 

9« 


tA*A»CtC> M0 ■ 199* Mr 

• 9»M^ • X 


00100 


c 




00IO< 

0* 


0l«t09tO0 



00lO« 

A0 


ifAklltl 



0OlO« 

7* 

( 




00100 

0« 


C0MM00 /9rkO0T/ 



0CIO0 

9* 



•0HfTOi|M0| 

.leoTAOlMPI .T|ST«f|or 1 

00100 

10* 


tcriari «c0(Mri 

•Cfirri 


00JO0 

II* 

c 




00100 

It* 


0Af« IS / 1 / 



00100 

19* 

c 




00100 

M* 

c 




00100 

10* 

c 


m 0Micw fi 

ttcltllp AktlTirOC kt(f 

00100 

I** 

c 




00100 

It* 

c 




00110 

!•* 


lri 0 *kt* 0 ktTAill>UO T 0 

000 


O0III 

10* 



001 


001 10 

90* 


1 • If 



00110 

91* 


trii.kt« 0 i.n 00 iti *00 to 

lOO 


001 It 

91* 


10 1 * 1 - 1 



00190 

99* 


lrt0«kV»0ktTA0lllt0O TO 

10 


00 lit 

90* 


kO TO 900 



OOltJ 

90* 


100 1*1*1 



0019* 



iriv0I*ALlfA0tlU*O TO 

IPO 


oei9» 

9t* 


1 • 1 • 1 



0Ol9t 

9I« 


900 If • 1 
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Ul • V • AtTT*M c I > 


00130 

3f4 

00130 

J0» 

00130 

314 

00130 

334 

00131 

334 

00137 

344 

00133 

3S4 

001 33 

344 

00133 

374 

00133 

3S4 

0013R 

344 

00I3R 

SO* 

0013R 

41# 

01 3R 

43# 

OOi3R 

434 

00I3S 

44* 

00134 

4S# 

00134 

44* 

00134 

47* 

00134 

4S4 

00434 

44* 

0013T 

SO* 

OOIMO 

SI* 

OOINI 

S3* 

OOMI 

S3* 

oom 

S4* 

00 Mi 

ss* 

004RI 

S44 

00M3 

S7» 

00 m3 

SB* 

OOMR 

S4# 

OOMS 

40* 

00M4 

41 4 

oomt 

43* 

OOlSO 

43* 

001 Si 

444 

00IS3 

4S* 

00IS3 

44* 

0CIS3 

47# 

OOIS3 

4B* 

00IS3 

44* 

00IS3 

70* 

OOlSR 

71* 

OOlSS 

734 

00IS4 

73* 

OOlST 

74» 

OQUO 

7S* 

00 U 1 

7*» 

00143 

77a 

00143 

7i* 

00143 

79* 

00163 

6 '• 

00I6R 

S M 


C 

C 

c 


c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 


CmPcOT CW^IC IHTEKPOtATlON TO COMPUTE P«CSSUPC/itA iCVEt PaT|0 

FanSUI ■ PPSTAOU’ ♦ lKP(3«t>«OX * CP < 3* I • 1 > > *0 X • cPf3«|-2M*0X 
FAnSCS) ■ FANS i I > /PPSTab I I » 

tMPLOy INtEPPOLATlON TO cOMPu^C 0EN5lTy,SCA Ltv£^ PaT?0 

FaNS(2) * Pn0Tab:|I ♦ HCP|3«I)*0X ♦ c" <3*1 *1 I I «DX « CPf3«r«2M«0X 
FaNSCAI ■ p %«IS f 2 ) /PHOTAB ( 1 I 

employ cu^tc INTCPPOLATiOn to compute speed of SOuND/SEA level 
Patio 

FANSI3} ■ SOMTABd) « (ICS(3*|I«0X ♦ cSf3«i*in*0X « cS(3»l*2n«0X 
FanS(T) ■ fans I 3>/S0MTAB I 1 ) 

employ CUbU INTEPPOLATiON to compute coefficient of vISCOSITT/SEa 
LEVEL RATIO 

Fans(R) • VISTABCn « UCy<3«t)«0X ♦ cv(3*|«l | )»0x ♦ CVC3*I-7M*0X 

FANSC8) ■ FANSIM) /V ISTA8 ( II 

PETuPN 


initial BqOnOaRY values - FUVdlO/SEA LE^EL PaT|0 


Sqo fans! I » 

FANSf2> 
fansoi 
Fans M I 
FanSISI 
fansui 
FANSITI 
FANSISI 
NETUPN 


• PRSTABI 1 I 

• RHOTAbUI 

• SONTABI I I 

• VtSTAbl I I 

• 1*0 

• I *0 

• I «0 
■ 1*0 


Terminal boundary values - function/sca level Ratio 
S ot fans ( I ) • prstab InPI 

FanSI 2> • RHOTABInPI 
FANSC 3 I • SONTABInP) 

FANSM» • VISTABINPI 
FaNSISI • FanSI I I/PRSTaBI I I 
fansui • FANSI 2I/RH0TA8U ) 

FANSI7 I • FanSISI/SOnTABI I I 
FANSISI • Fans ( A I / V I stab 1 1 1 


c 

c 


return 

ENO 


I 


tNU OF compilation: no diagnostics. 
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|l«« r««T««b * ^aiC It <»• -Idle* kCffCL |I|P1««I9A> 

»h% •«*( t#« II 7« 4T JP|W;I( 


sw««pvl|ht t^wm t«liT pplaT OVOPt* 


t'a***!. wPiPl CPPMIl 0007S7I DAttlOl OyOllOI OOCQCO 


»CI 


iIlBC** "IWII 


PPOi 

Mya fttkB)* 


»CMtON, BftnCI 


OMI 

9991 9» 

IL 

0901 

0901)9 

1 IM 

9001 

990 lai 

l)«« 

•901 

aoo)*i 

!*)• 

0001 

00n*7) 

l«IO 

•091 

99991' 

Ift’k 

9901 

O90d« 

1 

9901 

99f*lt< 

It 

•991 

fOOXO* 


90«l 

ooe«*i* 

foit 

0991 

•99191 


9901 

099 >9^ 


9001 

990)t« 

<0f 

•991 

<*909*9 

0*»L 

0001 

00'*l)l 

»0I. 

9991 

POO»*a 

'iL 

••90 ] 

1 00990 t 

1 

9900 

OOoOia 

|yj9t 

9009 

t '*90909 

m 

0000 ; 

1 O0000> 

ny 


99101 

1* 


ftiioto 

utiiit 

•9191 

t* 

c 



99l9> 

)• 


91*10)10* ll|ltVllf.8(tl«*ll>«CUl 

9919) 

a* 

c 



9019) 

%• 

c 



9919) 


c 

H 


•910) 

)• 

( 


■tPaiSC*! Tm( ruMrTiQv |m TnC OCtlOy Of Ocri«tT|Oa 

09l0< 

• • 

c 



9910) 

t* 

c 

« 


0910) 

10* 

c 


tALUll or f*c l*OC^t* 0 C«T Va«|*|LC I* ISCC*0I*0 0 * 9 f« 

•910) 

11* 

c 



0010) 

It* 

c 

f 

A Oai>oi*t>vieaik AiKftT ei*i*Slo*ce i* * cOaTatmiio Tr*c 

0010) 

O* 

c 


fAtwtl or Tm( ttCrt^DCK* ai«|*0LC cOO*t )rO»0 1 aa 1* 

0910) 

M* 

c 


icoucacc 10 I* twC ) «a0ta 

•910) 

10* 

c 



0910) 

la* 

c 

Jl .)* 


9010) 

17* 

c 



9919) 

10 * 

c 

Jl • 

1 • ri*)T Dia|vtT|vt gr T*<c runcTie* at TMf riaiT raouciiCB 

0910) 

la* 

c 


data roiaT ti t'vrrLlco i< Din 

eolol 

to* 

c 



00 10 ) 

ti* 

c 

yi • 

t - UcOhO 0C9tfAT|V[ or TyC ruacTfO* at Ih( rt«lt ToouLaTCD 

0010) 

at* 

c 


eart rolaT |i surrtico in OMi 

OOIM) 

t)* 

( 



9010) 

ta* 

c 

•il • 

) • Fiasr aab HC9N0 Bc«iaaT|rcS Or Tyc pumctiom at tmc ri0<T 

0910) 

»• 

c 


DAT* rofyi aAK furpLIIo 1* Oill a*0 Dttl. oCSrCCT 1 <rCLT 

•OlO) 

t»* 

c 



•Oiu) 

tr« 

c 

VI • 

a - M»st «MS aitOyO gt*l*|Tl*C) Or th[ rua-rio* at tmc u*if 

OOlO) 

t«* 

c 


gara roiNt a*l »wr»Lll8 1 •• Oili amo Oitl, *l *rcc T l att* 

00 lu) 

tt* 

c 



0010) 

)0* 

c 

Vn • 

1 - iiA)T DC*l>aiivi or (wiMCiina a' rnf tati iatuk,niD 

00 lu) 

il • 

c 


gala rgliil tl )urrLtCO 1 >. 01)i 
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t 


1 


oolua . 

32* 

c 



oolo:' 

33* 

c 


JH • 2 • StCOMO OEKivaTiVt or TmE ruwc7|0n *T TmC last TAaotaTtO 

CO t b 3 

i*« 

c 


data point 1$ supplied bn 0121 

00103 

3># 

c 



00103 

34« 

c 


0 • A ONt-OlHENSlONAt AKKaT OIHENSIoNEO EY 2 cOHTA|N|m6 THE 

00103 

37» 

• 


vaLoes or The oC»uvative 5 or the ruNciiON sPECirtro in 

00103 

33* 

c 


ACCoHOANCE NITn J2 and jn 

00103 

3f • 

c 


C • A On£-0|NENS|OnaL APKaY 0|HEN$|0NED 8Y 3#fH.n containing 

00103 

^0* 

c 


2hE C^sU 3^1-I**E I^TtProtA2|N6 1 c 1 t^TS fO»1 EaCM 

00103 

‘M* 

c 


interval 

00103 

12* 

c 



0QI03 

N3» 

c 


• - A OnC-OINENsiONAL array OIReNSIONEO «v n used as a RCRKING 

GOlOl 


c 


arrat only 

00103 


c 



00103 


c 



00103 


c 


copruTc The initial values or tjc recursive relations in 

OOlOJ 


c 


accordance rith The SPcirieo proRcRTy of the ruNctioN at the 

00103 


c 


initial BoUNOARY 

00103 

&0* 

c 



00103 

Sf 

c 



00)0*1 

hz^ 



so to 1 1 •2»3.R I • J| 

OOlOS 

S3» 


1 

• in ■ -o^s 

0010^ 

Si* 



cm • 3»0/M(2) - jm>l*Hr<2l - Yim/IxlEl • AfiM - Ofin 

00103 

%%• 



so To 20 

001 10 

%k^ 


2 

Rin ■ Oto ^ 

001 11 

S7* 



cm • om 

00113 

sa* 



so to 20 

00113 

sa* 


3 

Rill • 1*0 

001 M 

ao* 



cm • 3*o*o<2) * a«o/iiC 2 ) • xim*(D(ii * <ri 2 i • rmi/iK( 2 i * 

OOU^ 

ai* 



00I1& 

a2* 


20 

N • N - 1 

0011^ 

a3* 

c 



00M» 

aa* 

c 


coHPurc all inTerhioiaTc recursive relation values 

oons 

as* 

c 



0011* 

aa* 



00 30 1 • 2.H 

00121 

hf 



• m • »ixii*i) • xim/( 2 i 0 *fx<i*i 1 • xci*m ♦ ixm • xm-ii.*» 

00121 

aa* 



xii^m 

00122 

at* 


30 

cm • Ri|)/(xci*ii • X m 1 «f a*o*i 1 vm - ni-m/ixm • xii-in 

00122 

ro* 



A- iYii*i) - Ymi/(i(i«ii • xiim ♦ ix(ii • x(t»iM*c(t«m 

00122 

2|* 

c 



00122 

72* 

c 


coHpuTt The terminal values or -he Recursive relation in 

00122 

2J* 

c 


accordance RITm The SPLClrlEO property or the rUNcTlON aT the 

00122 

Ti* 

c 


Terminal boundary 

00122 

7S* 

c 



0012M 

7a* 



Ir 1 j) •sT* 2 IG 0 To 70 

0012* 

7;* 



IF ( JN,E0*21S0 TO *0 

00130 

7a* 



CINI • (6,0/(I(NI - ||N«tn*(0(2) - lYiNl * Y|N«m/(AlNl « XlN*H 

00130 

7t • 



aM « ClM))/(2«0 * R(Nl) 

00130 

ao* 

c 



00131 

• i* 



sc ro so 

00132 

i2* 


*0 

C(H) « 0(21 

00132 

83* 

c 



00132 

at* 

c 


COMPUTE The aRffoximaTe second derivative Values or Tme function 

00132 

as* 

c 


AT the srtCIFIEO DATA POINTS 

Q0i32 

a** 

c 



COI 33 

8M 



GO TO SO 

OOli** 

a«« 


70 

UO 7 2 1 « 1 k N 

00137 

at* 


72 

C t 1 *N 1 • C ( 1 1 


ORIGINAL TAv... 
OF POOR QbAA.I] 
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cn*MM/«i 1 1 


OOI^I 

*0* 


Cl 1 1 • 0121 

ooi^a 

♦I* 


00 33 1 • liH 

oon» 

f2* 


73 C( 1*1 1 • ICC 1 1 . Cl |/fl| II 

00l«i7 

tj* 


40 TO 31 

00 ISO 



«l n « N - 1 

OQISI 



•|M| ■ 1*0 

00IS2 

»4* 


CIN) • 0121 

00IS3 



CM! • 4*0/UlNl * 2|Hn*<Ulll * CT(N| • TIN))/|X|N| - XIhIII - 

OOlSi 

to* 


X3«0»oi 21 

OOlSH 

• 


HH • H • 1 

OOlSS 

100* 


IF INNI |S0« 

001*0 

101* 


00 to 1 * 

001*3 

102* 


Mill • -l•0/^xll•|l - Emi*M*ll*2l • K« 1* 1 M/«< 1 ♦! 1 ♦ 2.0«lxl!*t 

001*3 

103« 


- 111 Ml 

001*^ 

I0*« 


•0 Clll ■ I•0/IIC|*1I - XC 1 1 l•(C( 1*1 l*lll|*2l « XI |»i ! !/■! 1*1 1 - 4.0* 

00l**» 

US* 


iiiYii*ii • Tini/ixi|*n • xiiii • iYii«2i • Y( i«i n/ixf 1*21 - xn 

001*^ 

104* 


X»1 1 1 1 1 

OOi** 

IU3« 


so DO *0 1 • rttl»*l 

0013 1 

loa* 


40 Clll • c«ii ♦ •ui*:ii*ii 

00133 

lot* 


71 00 41 1 • 1 .N 

0013* 

no* 


41 *111 • c< 1 1 

0013* 

III* 

C 


0017* 

112* 

C 

CONFUTE ThC CUBIC spline I k rCBPOLxT t OH cOerrtClENTS rOR The I«STi 

0013* 

113* 

c 

2«N0« and 3*110 OeSNEt TCRHSt IN ThAT OBOEB 

0013* 

II** 

c 


002U0 

IIS* 


DO 42 1 • ItN 

00203 

114* 


CI3*|»2I • Irimi * T1 1 1 1/ IXI I *t 1 • X(|ll • l■l|«l) ♦ 2fO*BCtl’*l 

0Q203 

113* 


|X| 1*1 1 • III 1 1/*«0 « 

0020* 

1 in* 


C|3»|«l 1 • 0»s«4l 1 I* 

0020* 

II** 


42 ci3«i) • iBini) * 01 1 1 )/i6«o*ixi i«i 1 • xinn 

0020* 

120* 

c 


00203 

I2i* 


RktuBn 

00210 

122* 


Eno 


tNO 0^ CONMlATIONI mo 0|*6n0STICS. 
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v«i|*Ai tlO» rb«r«f t II kfvtw f** it|0|Ooio»i 
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OOlOl 
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kPI 1 • 1 , NfM • 1 If , 
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eolo' 
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0010’ 
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0(21 
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u* 
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c 
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*Cf t**> 



001 1 1 

1 7* 
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1 4* 


D4 f4 tALfkSfl 1 • . I***! 1 
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c 
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i 

f 4 
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7>4t2. t>0tl* 

1 .1VC1> 

1 

001 11 
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c 

7 

1 

f 

13 11 

11 


001 11 
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i-’vei. 

1.0li> 

l>14tj* 2*4Cl* 

iwsci. 


001 ll 

ll“ 

c 

11 

1 • 

14 

1* 1 ’ 

14 


00 111 

If* 


M>0il. 

1»*C 1. 

f.Of 1> 

«>4ti. 4-OCl* 

4.011. 


001 |1 

li“ 

c 

I* 

10 

ll 

12 11 

If 


001 ll 

l4“ 


If.oil. 

i*oii* 

*.0l 1 > 

1 >01* . 1 • l(f » 

1 .ICf * 


001 ll 

17* 

c 

t% 

2* 

27 

!• 2« 

10 


001 |1 

14* 


M . n*i 

1 • f t f • 

1 ,4Cf • 

l»4lf, l-IC** 

1 .11* > 


001 |l 

If* 

€ 

1 1 

J1 

li 

If IS 

34 


001 1 1 

In* 


1 1 • «i<* . 

" ^ ^ 0 

l.ltf • 

/••iff 2»*|f* 

2*0l*. 


001 |1 
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